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ABSTRACT: We calculate the triviality bound on the Higgs mass in scalar field theory 
models whose global symmetry group SU{2)i x 5'C/(2)(,yg|;ojijg^i ^ 0{4) has been replaced by 
0{N) and has been taken to infinity. Limits on observable cutoff effects at four percent 
in several regularized models with tunable couplings in the bare action yield triviality 
bounds displaying a large degree of universality. Extrapolating from = oo to = 4 
we conservatively estimate that a Higgs particle with mass up to 0.750 TeV and width 
up to 0.290 TeV is realizable without large cutoff effects, indicating that strong scalar self 
interactions in the standard model are not ruled out. 



1. INTRODUCTION. 



This paper examines the regularization scheme dependence of the triviahty bound on 
the Higgs mass in a 0{N) symmetric scalar field theory to leading order in 1/N. Our 
purpose in doing this is two- fold: we wish to investigate the issue by explicit, analytical 
calculations and we need the results both directly and indirectly to complement numerical 
work at = 4. The results are needed directly for estimating cutoff effects on physical 
observables that are not accessible by Monte Carlo and indirectly for guiding our search 
in the space of lattice actions to the region where heavier Higgs particles are possible. 

A preliminary account of some of our results has been presented in [1,2,3]. Our general 
reasons for suspecting that the present numbers for the bound are too low in the context 
of the minimal standard model have been explained before [1,4,5] and will not be repeated 
here. 

The basic logic of our approach has also been explained before [6,5,7] but, in order 
to make this paper more or less self contained, we shall briefly review the main points. 
We are working with a scalar fleld theory at a finite cutoff. The field has N components 
(A^ = 4 for the minimal standard model) and the action is 0{N) symmetric. We are 
interested in the broken phase in which the symmetry is spontaneously broken and wish 
to pick cutoff schemes that preserve as many of the continuum space-time symmetries 
as possible. Our cutoff scheme must also exclude unitarity violations at the energies and 
momenta of interest which are small relative to the cutoff. 

The physical scale is set by Fyr, the "pion decay constant", and we are interested in 
Mjj the "Higgs mass'' . Because of triviality the ultraviolet cutoff cannot be removed and 
therefore cutoff effects are unavoidable. We require those cutoff effects that are observable 
to be bounded by a few percent in relative magnitude. This imposes a limit on the ratio 
Mjj/Ft^. Generically, the leading cutoff effects are of order inverse cutoff squared. It is 
possible in some models to arrange that the leading cutoff effects be of order inverse cutoff 
to the fourth power instead. We regard such "finely tuned" models as overly contrived. 
The cutoff effects of the low energy effective theory are ultimately determined by the 
underlying full theory and we believe it is highly unlikely that such a "finely tuned" model 
would be produced without any particular reason. We wish to find a region in the space 
of actions where, without excessive "fine tuning",* Mjj/Ft^ gets to be as large as possible. 

Since the cutoff effects are assumed to be small and "fine tuning" has been excluded we 
conclude that we can use the RG to count the dimensions of the space of actions we ought 
to look at. We have one relevant operator of dimension 2, two operators of dimension 4 
and four operators of dimension 6. We do not need to consider higher dimension operators 

* Except the one possibly needed to make F^^ small relative to the cutoff in the first place. 
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because their effects are subleading and we have excluded the fine tuning of the leading 
effects to zero. Since wc are interested only in physical effects, namely in cutoff effects in 
the S-matrix, we have to eliminate the redundant operators induced by field redefinitions. 
For dimensions up to six there are three such operators. Thus, we end up with a 7- 
3=4 dimensional space. One dimension corresponds to the relevant direction and sets the 
scale, another to the marginally irrelevant parameter and the two remaining parameters 
span all the possible leading order observable cutoff effects. It is preferable to think only 
about dimensionless physical observables with the natural physical quantity obtained by 
multiplication by the appropriate power of F-^. All momenta are also measured in units of 
and so is Mjj; the marginally relevant direction can be thought of as parametrized by 
g = ZMl/Fl 

Intuitively, we would expect to increase g by approaching the nonlinear limit where 

— * 

the field has fixed length. Hence, for mass bound studies one can restrict one's attention 
to nonlinear models. This is known to be true with naive lattice actions and we have 
checked it for some other actions N = oo. Hence we shall concentrate on nonlinear 
actions. At the most naive level one can think of a nonlinear action as obtained from 
a linear one by taking the bare four point coupling to infinity and adjusting the other 
couplings so that the model has a limit. From this point of view one of the four free 
parameters has been eliminated; in other words, the effort to maximize g is assumed to 
guarantee the elimination of one of the four parameters, leaving us with only three. The 
nonlinear actions can also be viewed as chiral effective actions in the sense defined by 
Weinberg [8] . The low momentum behavior of any broken theory (even with a cutoff if the 
latter is sufficiently symmetric) can be parametrized by such an effective action and, to 
first nontrivial order, one needs three parameters. This is exactly the number we came up 
with after reducing by one the original dimensionality. Therefore we know that the three 
parameter set of nonlinear actions we intend to focus on is reasonably general, since it can 
at least reproduce all "pion" scatterings to first subleading order in the external momenta. 
To be sure, we still expect the cutoff effects to depend on two parameters as before the 
restriction to nonlinear actions, because such a constraint should not reduce the number 
of contributing operators in the vicinity of the Gaussian fixed point. 

Our further approximation will be to take N to infinity. We shall see that this effec- 
tively reduces the number of parameters from three to two and, moreover, at = oo the 
cutoff effects we shall be interested in can now be parametrized by a single parameter.* 
The situation becomes as simple as it could get: we have one parameter that sets the 
scale ( "coarseness" of grid in the lattice case) and another that represents the entire free- 
dom that is available at first subleading order in the expansion of the S-matrix in inverse 

* Even a bare action that has only a single free parameter, will, in principle, have observable 
cutoff effects that are parameterizable by two coefficients, associated with the leading 
"irrelevant" operators after elimination of "redundants" . 
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powers of the cutoff. For fixed value of the second parameter one can imagine that some 
extremization of g has aheady been carried out by taking the bare four point couphng 
to infinity. This simphfication at infinite N makes it easy to relate very different cutoff 
schemes and leads to a reasonably "universal" bound on g. These facts, derived below, 
provide a basis confirming the general validity of the N — oo results^ obtained previously 
by Einhorn [9,10]. 

We shall work with a class of Pauli-Villars regularizations parametrized in addition to 
the continuous parameters discussed above by an integer n > 3, and with lattice regulariza- 
tions on hypercubic and F4 lattices. When the range of the additional couplings is suitably 
restricted (but not "finely tuned" ) all these regularizations give very consistent estimates 
for the N — 00 bound. We present our results properly scaled to = 4 and obtain, with 
Ft, = 0.246 TeV, < 0.820 TeV with rather stringent bounds on the cutoff effects; 
changing the latter by a factor of 10 may affect the bound by about 0.050 TeV in the obvi- 
ous direction. We can try to guess what correction on the bound might be due to N being 
four rather than infinity by looking at the available numerical data at AT = 4. This leads 
us to expect for A?" = 4 a bound somewhere between 0.750 TeV and 0.800 TeV. While 
the relatively recent previous estimates were not much lower* (0.600 TeV to 0.650 TeV) 
the effect on the width, F^, is more significant: Fjj/Mjj may reach 0.4 when the bound 
is close to saturation and this implies that strong scalar self-interactions without strong 
observable cutoff effects are possible. 

Since the previous bounds were obtained from a very restricted class of actions, chosen 
just because they were simpler to analyze for technical reasons, it was somewhat premature 
to put forward the 0.550 TeV to 0.650 TeV range as the lattice triviality bounds with 
lattice spacings between l/{5Mjj) and 1/(2M^) and imply direct relevance to experiment 
[13]. To be sure, we do not contest the validity of the numbers within the particular 
regularization that they were obtained in, and, in retrospect, they were not very far off even 
in the general context. Our analysis shows that we are nowhere near a ratio Mjj/F-^- ~ 6 
that would be expected in a QCD like theory and our Higgs is quite "elementary". We 
have no explanation for why the range 3 < Mjj/F-,, < 5 seems to be so difficult to attain 
with "an elementary" Higgs; this looks like a stronger numerical reflection of triviality 
than one might have originally suspected. 

Throughout the paper we shall use a notation with capital letters, Ftt, Mjj and Fj;/, to 
denote the physical quantities, measured in TeF, and with lower case letters, /t^, mjj and 
7^, to denote those quantities in units of the cutoff. A, or for lattice regularizations. 

t With some small, inessential corrections. 

* Very early estimates [11,12] were actually around 0.800 TeV but this was accidental; these 
authors also used the simplest possible lattice action for which we have nowadays better 
results. 
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In the large N limit the renormalized vacuum expectation value (pion decay constant), /tt, 
diverges as a/ZV but and 7^ stay finite. We therefore rescale fj^ and F-j^ by ^/N to 
make them also finite at — > cxo. From now on and until sections 7 and 8 /tt and F^^ refer 
to the rescaled values. Of course, when we use the large results for = 4, we undo 
this rescaling; when we present our large N results in physical units we take the rescaled 
Ftt as 0.246/^4 TeV = 0.123 TeV. 

The outline of the paper is as follows: In the next section we shall briefiy analyze 
linear models with particular higher dimensional operators to see explicitly that in order 
to estimate the bound we have to look only at the nonlinear limit. We proceed to define the 
class of nonlinear models we shall consider in the rest of the paper. In section 3 we work 
out the large phase diagram of the models in the Pauli-Villars class. Section 4 deals 
with the cutoff effects in the Pauli-Villars class of models, first in a simplified manner and 
then in detail. Sections 5 and 6 repeat the analysis for the lattice models. In section 7 we 
present the results of our large N study of cutoff effects for the Higgs width to mass ratio 
and for Goldstone pion scattering and the implications on the Higgs mass bound. Also, a 
comparison between large N and available numerical results at A?" = 4 is carried out. This 
allows us to make inferences from the results obtained in this paper to the AT = 4 case 
of actions that have not yet been studied numerically. Our conclusions are presented in 
section 8. Three appendices have been added to present some technical details and expand 
on side issues. 



2. LINEAR VERSUS NONLINEAR ACTION. 



In this section we briefiy discuss the linear case. We shall see why we were led to 
consider only nonlinear actions in the sequel. Also, by comparing the formulae derived 
later on for the nonlinear case to the ones in this section, one can see explicitly that the 
nonlinearity of the bare action has no observable effect in the critical regime as expected 
by general RG arguments. More precisely, by making physical measurements in the critical 
regime and at low enough momenta one cannot say whether the bare theory was linear or 
not. It is because of this that simulations employing a nonlinear action have something to 
say about the standard model as we know it. 

We shall use the following notations: The partition function is Z, the metric is Eu- 
clidean and the action is denoted by S. Z is given by 

Z = J exp(-5) , (2.1) 

where the fields have been suppressed in the action and integration measure. For spatial 
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{x,y ■ ■ •) and momentum (p,q,k---) integration we use 



and as usual 

A denotes the cutoff. We shall suppress arguments whenever possible. 



(2.2) 



(2.3) 



2.1. The general linear action and the saddle point equations. 



According to the discussion in the introduction the most natural set of actions we 
should be studying is given by 



S 



/ 



(2,1.1) 

The subscript "0" denotes dimensionful parameters. The kinetic term is regulated in a yet 
unspecified manner, but it is assumed that all subsequent potential ultraviolet divergences 
are cut off at about A. Thus we write 



(2.1.2) 



(2.1.2) contains seven free parameters (two of them are implicit in Ko(~^o))' eoch. 
operator of dimension < 6. Later we shall choose to fix some of these parameters to 
simplify our analysis. 



We scale A out of the problem by defining 



X 



= Aa;o , 00 = A0 , /^o = ^^/i^ , mfi = ^ , ^2,0 = "p , ^3,0 = ^ (2.1.3) 

(2.1.4) 



m 



s = 



and all our variables become dimensionless 
^1^ ..2 



4N 



2N 



2N 



Introducing into the functional integral 



da{x)— / dw{x)exp 
-oo 4^ J—ioo 



N 
6 



u>{x) 



a{x) - 



(2.1.5) 



we obtain a new action with more fields 



Si = ]-j$[K- md^ad^ + u;]$+^ f [r]2{d^,af + i/a + + - au] . (2.1.6) 

— * 

The zero mode of needs special treatment requiring that we separate it out explicitly 

$(x) = ^/Nv + vH{x) + tt{x) , '0=-, \v\ = V , ■0-7r = 0, / H{x) = / 7t-'{x) = . 

^ Jx Jx 

(2.1.7) 

We integrate out the H and tt fields and, with K = K — ijid/jadu + a;, restricted to 
operate in the space of non-constant functions, we obtain 



J J X J X 

(2.1.8) 

Here uj' is the non-constant part of u 

d^^{u - cj') = , [uj' = 0. (2.1.9) 
Jx 

u)' rather than u comes in because coH = co'H as a result of the zero mode, v being 
separated out. In the functional integral one is still left with a piece v^~^dv because 
V hasn't been integrated over. In the infinite volume limit, and assuming translational 
invariance, this piece can be ignored.* Taking N to infinity we obtain the general saddle 
point equations 

f 1 2 

/ T^r 9n 9 1" ^ = 

2 (^.l.lUj 
/ P , 2 , X , 2 

We notice that 772 has disappeared from the saddle point equations because we assumed 
translational invariance of the saddle. 

If we keep w as a free variable and plug in the solutions of the saddle point equations 
into S2 we obtain the effective potential. Here, however, we just wish to fix v at the 

— * 

vacuum expectation value of 0. In the broken phase one then has a;s = 0. For > ^ v 
will vanish. 



* The large N analysis can be easily extended to finite volumes where the probability distri- 
bution of V can be calculated showing explicitly how the infinite volume limit is approached. 
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We are looking for the candidate transition surface. It is only "candidate" because 
it might be cut by other transition surfaces rendering some portions of it metastable. 
However, any true second order transition point we are interested in must be on this 
surface. The condition for criticality fixes ji^ = iJ>c{r]i, A). //^ is implicitly defined by 



I 



1 

~ o - (^s 



+ ^^^^^^ 

For the time being we ignore questions regarding global stability. We take 5/1 = /i —/i^ 
to be negative and small relative to unity. This should place us in the broken phase close 
to the transition. 



2.2. Small fluctuations around the saddle point in the broken phase. 



We now wish to compute the propagators. We take v to point in direction 1 in the 
internal space. Hence tt^ = 0; we use latin letters to label the N — 1 nonvanishing com- 
ponents and view the "pion" field as made out of them only. The N = oo propagator is 
immediately read off, in Fourier space, as 

xab xab 
<7rV>=— ^-^ ? = • (2-2.1) 

To get the Higgs field propagator we integrate out the fields 

^3 = ^^-^TrlogK + ^ [ HKH + [ u + VNv I uH+ 



^ /"r /CI \2 , 2 A 2 , % 3 1 

y J midfxCT) +iia + -a + —a - auj] . 



(2.2.2) 



At large N we expand around the saddle point with Sa = a — as, Suj = uj and SH = H. 
We get, neglecting order 1/N corrections, 

= ^ [ [l-5aBaaSo- + hujB^ujSu; + 6aBau;Su;] + ]- j SHAgSH + VNv f 5u5H+ 

^ Jp 2, Z Z Jp Jp 

(2.2.3) 



where 

As(^p-Q)A5(^p + g) 

To find the propagator of H one has to invert a 3 x 3 matrix whose entries can be read off 
from (2.2.3, 2.2.4). Note that 772 reappeared and the Higgs mass wiii depend on it. 



2.3. Restriction to 771 = 773 = 0. 



The problem simplifies considerably when we set ?7i = 773 = because 772 doesn't 
appear in the saddle point equations and affects only the Higgs mass. Note that using 
field redefinitions and keeping terms up to order inverse cutoff square 771 and 773 can be 
eliminated from the action (2.1.4). Intuitively 772 looks like a parameter that may have a 
significant effect because in (2.1.4) it gives some extra "stiffness" to the modulus of (p. We 
wish to sec whether keeping 772 > can lead to a situation where the Higgs mass is no 
longer monotonically increasing with A (since cr is a real field [see (2.1.5)] 772 must now be 
non-negative to keep the action bounded from below). We shall see that any reasonable 
limitation on the cutoff effects in tt tt scattering prohibits this from happening. Hence, 
as far as the triviality bound is concerned, we end up being driven to A = 00 and the 
dependence on 7/2 ultimately drops out from the Higgs mass too. We shall assume later 
on that turning on the other 7/ couplings would have had similar consequences, and that 
the triviality bound would be independent of them too. 

Our action in dimensionless variable is now given by 



S = 



(2.3.1) 



The saddle point equations become 

+ / r^/ 2! , = (^s , + AcTs = Us . (2.3.2) 

Jp K{p^) + UJs 

One can show that there are no competing saddles here and hence, once 771 = 773 = 0, 
questions about global stability do not arise. 

The critical line is given by 



1 



and the broken phase is where 



There, u>s = and the vacuum expectation value v is given by 



(2.3.4) 



(2.3.5) 



Denoting the fields SH, Su and Sa by 1^2 and il^^ respectively, we obtain the 
following propagators in Fourier space 



where 



M 







N 



N 

4(A + 2r72p2)y 



and 



laKav-q)Kih 



(2.3.6) 



(2.3.7) 



(2.3.8) 



Instead of looking at the Higgs propagator we look at tt tt scattering. The invariant 
amplitude M for the process 7r''(l) + tv^{2) 7r^(3) + 7r^(4) is 

out < cd\M\ab >in = A{s, t, u)SabScd + ^t, s, u)6achd + M^, t, s)6adSbc 



s = {pi+P2) , t = {pi-P2y 



u 



{pi - pa) , P 







(2.3.9) 



in standard Minkowski space notation. At N = 00, A{s,t,u) depends only on s and, in 
our restricted model, is proportional to the Suj propagator at p'^ = —s. Hence, by looking 
at the Slij propagator we can get both the Higgs mass and estimate the cutoff effects in tt 
TT scattering and on the width to mass ratio for the Higgs. From (2.3.6, 2.3.7) we obtain 

2 K(p2) 



(2.3.10) 



We now make our choice for K{p^): 

K{p^ 



P^ + ip^f + <pY 



(2.3.11) 



By setting the coefficient of the (p^)^ term to one we have restricted another parameter, 
leaving us with 3 free parameters. With a rescaling of the cutoff and the fields* one can 



* The rescaling of the fields, by itself, would lead to unphysical cutoff effects only, but the 
rescaling of the cutoff does induce a physically observable change in the action. 
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easily adapt the analysis to the case where the coefficient of the (p^)^ term is different from 
one, as long as it stays positive. We don't expect small changes in this coefficient to affect 
our final conclusion in any way, but have not checked explicitly whether our conclusion 
still holds when the changes are allowed to become large. 

To make //^ finite we need e > 0. We introduce an e dependence in //^ such that in 
the limit e ^ /j? diverges and Sfj.'^ stays finite. Everywhere else only 6fi^ appears and 
e can be set to zero.^ Thus, in the definition of I{p^) we take 



(2.3.12) 



which gives us 



lip' 



87r2 



-- log(p2) + (1 + ^) log(l _ 



1 + ^ arcsinh | ^ ^ 



(2.3.13) 



It is easy to check that the local stability requirement /(p^) > of the saddle holds for all 
Euclidean momenta; it is known [14] that if one extrapolates the leading asymptotic ex- 
pression for I{p'^) valid for small momenta to large ones, the stability requirement appears 
to be violated; this gives rise to the "tachyon" problem in the renormalized expressions and 
signals that for high momenta one cannot remove the cutoff dependence from the theory. 

When we analytically continue to the physical regime for tt tt scattering we replace 
by —w — zO"*" with positive w. We obtain 



I{-w - iO+) 



IGtt 
1 



Stt^ 



11 4 

- log{w) + ( 1) log(l — w) + \ 1 arcsin 

2 w \ w 



(2.3.14) 



2.4. Higgs mass in the restricted model. 

As a measure of the Higgs mass we choose the more accessible quantity, m^, defined 
as the smallest positive root of Re[< dudu It makes physical sense to only allow 

center of mass energies for which ghosts cannot be produced. On the other hand, the range 
of center of mass energies must be allowed to surpass mji by a factor of 2-4. A pair of 

^ Note that when the coefficient of the (p^)^ term is negative e cannot be taken to zero; 
nearest neighbor lattice actions have a negative coefficient for (p^)^ and therefore lattice 
actions ought to be analyzed separately. An example of such an analysis is presented in 
Appendix C. 
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ghosts can be created with w > A. We therefore restrict mji by mji < 0.5. Later on we 
shall sharpen this requirement somewhat. 

It is clear that a zero will develop in < SuSu > for w = A/^, where we introduced 
2^72 = C > (see eq. (2.3.10)). This is a cutoff effect of large relative magnitude that 
has to be forbidden at least for energies as low as w < 1 and maybe up to = 4. The 
unwanted zero is pushed to sufficiently high energies if we impose the requirement A/^ > 1 
or maybe even A/^ > 4. We shall be able to derive our main conclusion even with the 
milder restriction A/^ > 1 and therefore we shall stick to it from now on. 

Defining x = m|j and G = m^/(2v^) we have to solve the equation 



Stt 



2 



G 



= {l-x) 



Stt 



2 



X-^x 



ct>{x) 



4>'{x) 



X 



(2.4.1) 



with 



1 1 — X A f x 

(j){x) = - log{x) + —— log(l - x) + W - - 1 arcsin ( W - 

1 , 5 19 2 23 Q 251 4 

~ — logx H X H X H X H X + ■ ■ ■ . 

2 ^ 12 120 280 5040 



(2.4.2) 



Note the absence of a constant term in (2.4.2). When x is very small compared to unity 
we get 

Stt^ 87r2 1 



and clearly the maximal G is obtained at A = cxo. This conclusion has been obtained in 
the perturbative regime and applies only when G is quite small even at its largest value. 
We want to show that the conclusion holds irrespectively of the magnitude of G. 

Observe that (t>'{x) will have a maximum, < a;(A,^) < 1 and a solution x exists for 
a given v only if v satisfies (j) {x) > IGtt v . We are interested in the smallest positive 
root of (2.4.1); it satisfies a; < ^. A short analysis shows that x(A, ^) > a:(oo, 0) ^ 0.18 and 
that ^'{x)lx decreases monotonically for x e (0,C). C is some number that can be shown 
to be larger than x(oo, 0). 

The cutoff effects can be chosen to be characterized by m^. We now decide to limit the 
cutoff effects by the bound < y^p. We pick some positive p satisfying p < a;(cxD, 0) ~ 
0.18. We shall see later on that this restriction on p is not at all severe (we already argued 
above for a p < 0.25). The smaller p is chosen to be, the more stringently the cutoff effects 
are limited. 

One can view v^. A, and ^ as free positive parameters restricted only by A/^ > 1. 
Indeed, any desired value for can be attained by tuning ix^ in (2.3.5). We now look for 
the largest possible coupling G that satisfies (2.4.1) with some x < p. More explicitly, we 
are looking for the set {v^,X,^} for which the equation IBtt^v^ = (j)'{x) (see eq. (2.4.1)) 
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gives an x that leads to the largest G possible while x is not permitted to exceed p. Suppose 
we have found this set and its associated x and G. The monotonic decrease of (j)'{x)/x and 
eq. (2.4.1) imply 

^>^>-(l-p)^(p) . (2.4.4) 
G p 

For the largest G possible both inequalities in eq. (2.4.4) ought to become equalities. Note 
that (j) does not depend on A and ^. The first inequality becomes an equality simply by 
setting X = p, which shows explicitly that indeed the coupling is maximized when the 
cutoff effects are as large as they are allowed to get. This shows that, as expected, in order 
to make the coupling as large as possible one has to allow the cutoff effects to grow up 
to the bound one has declared from the beginning; in other words, the cutoff effects as 
measured by the mass to cutoff ratio and the coupling, defined by the mass to v ratio, are 
monotonically related. The second inequality in eq. (2.4.4) can be made into an equality 
only when A — > oo; ^ is not restricted as long as it is assumed to always obey < ^ < A. 
With any finite A > ^ > the coupling G would have been smaller for any x^ including 
the "best" value, namely x = p. Therefore the Higgs mass bound is obtained in the limit 
A — > oo and does not depend on ^ and hence the coupling 772. 

One can compute the true pole location and one gets, for example, with M^/A = 
^ = ^ = 0.397 Mh/A = 0.188 and T^/A = 0.198 showing that our restriction on p 
did allow for sufficiently wide (and hence heavy) Higgs particles. Note that Mjj/A and 

are very different numbers for such large couplings; still, in the whole range they are 
monotonically related so our main conclusion is unaffected. 



2.5. Nonlinear actions. 



In the introduction we counted the number of parameters our space of actions should 
depend on. In this section we saw that just having the right number of parameters is not 
sufficient; for example, with rji = rj^ = 0, varying 772 fhe physical range had no effect on 
the bound, which is obtained in the limit A 00. We also considered a linear model on 
an F4 lattice having the equivalent of the couplings 172 a-nd 773 and found again that the 
Higgs mass bound is obtained in the limit A — > cxd. This analysis is sketched in Appendix 
C, since it relies heavily on methods developed in subsequent sections. We believe that 
the same conclusion would be found if ?7i, 772 and 773 were allowed to vary simultaneously. 
When A — > 00 with a properly adjusted ^j? — > —00 we get a nonlinear action and the r/i 2,3 
either drop out or can be absorbed in K. In the nonlinear limit, power counting reduces 
to derivative counting. The leading cutoff effects are now parametrized by "dimension" 
four operators (the leading nonlinear term is of "dimension" two), and we therefore need 
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to add four-derivative terms to the non-linear action. There are three such terms, but one 
is redundant and can be ehminated by a suitable field redefinition. 

We now have to address the question whether this time there will be some dependence 
on the couplings when they vary in reasonable intervals. By "reasonable" we mean intervals 
of order unity for the dimensionless couplings; we expect all sorts of effects to limit the 
range in which these couplings can vary but it is almost certain that when everything 
is taken into account the allowed intervals will be of order unity. We now argue that 
indeed, with a nonlinear action, sizable variations can be induced in physical observables 
by variations of order one in the couplings. 

For this purpose it is useful to view the nonlinear action as a chiral effective Lagrangian 
rather than a cutoff version of the usual renormalizable scalar field theory. Let us do some 
numerology: From previous work we know that for the bound we have Mjj / {2Fty) ~ 2.5* 
and A, the maximal allowed "momentum", is about 2'kMu. Thus K/ {At{{2Ft^)) ~ 1.25 and 
loop effects will be of relative order one for observables depending on external momenta 
smaller than A [15].^ Since loop effects are more or less of the same order as tree level 
effects, we can reasonably expect that order one variations in the bare couplings will have 
measurable effects. In short, unlike in the linear case, investigating the dependence of the 
bound on the additional terms won't be a waste of time. 

All our actions have a derivative expansion of the form 

Sc = 

(2.5.1) 

where ^ = Nj3. The parameter can be absorbed in hi and 62 by a field redefinition 

- ^_tM!L_ . (2.5.2) 



^02 + 52(^20)2 _ 2600^20 



Out of the four parameters /9, 69, 61 and 62 oiily three combinations affect the leading and 
subleading terms in the expansion of pion scattering amplitudes in the external momenta. 
With Pauli-Villars cutoff we shall simply set 60 = 0, but on the lattice it is usually more 
convenient to stick with the particular value for 60 that comes from the expansion of the 
lattice kinetic energy term at small momenta. This has to be taken into account when 
a bare action regulated by Pauli-Villars propagators is compared to one regulated by a 
lattice. At infinite N the 62 term won't contribute to the saddle point equations because 



* Recall that we are working with a rescaled F-^. 

^ There are arguments that would lead us to consider, in our normalization. A/ {A'k{-\/2Ft^) 
1.76 instead [16]. 
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the 0{N) invariant bilinear that appears in it cannot acquire a vacuum expectation value 
as long as Lorentz invariance is preserved. More analysis will show that this term doesn't 
contribute to the leading cutoff effects on tt tt scattering and the Higgs width. Hence, up 
to questions of global stability 62 can be ignored and simply set to zero. This fact will 
be used as an argument to simplify the kind of lattice actions we are going to consider 
explicitly. 

The class of actions regulated by Pauli-Villars terms we choose to investigate is, in 
dimensionless units, given by 

S = 

(2.5.3) 

where = Nj3 and n > 3 is an integer. When n ^ 00 we approach a sharp momentum 
cutoff. For such a cutoff the parameter counting presented in the introduction doesn't work 
because even at first subleading order in the inverse cutoff nonlocal higher dimensional 
operators make their appearance. We are interested in large n values, however, because 
there one expects a greater similarity to the lattice, in the sense that arbitrarily high 
momentum excitations are almost totally suppressed. 

We shall also study actions regularized on the lattice because only the latter can 
be solved nonperturbatively by Monte Carlo. We shall consider the F4 lattice and the 
hypercubic lattice. On a hypercubic lattice the operator counting is slightly wrong because 
of the existence of a Lorentz breaking dimension six operator which has no counterpart 
in any conceivable extension of the minimal standard model. This has no effect on the 
particular obscrvables that we will consider at = 00, but would lead to measurable 
effects at order inverse cutoff square in other observables. Such a problem does not arise 
for the F4 lattice. 

The most concise action on the F4 lattice that reduces, up to higher order terms, in 
the continuum limit to the form (2.5.1) is 

S= -2Npo J2 ^{x) - Six') - Npi J2 [^ix)-^{x')f 
<x,x'> <x,x'> 

X <ll'> 

inxVC i'ns'VC, x,x',x" all n.n. 

(2.5.4) 

where x,x',x" denote sites and < x,x' >^l,l' links. Here the field is constrained by 
^"^{x) = 1. Geometrically, the first two terms couple two fields connected by a nearest 
neighbor bond and the last term couples three fields that live at the corners of an equilateral 
triangle whose sides are nearest neighbor bonds. The existence of elementary triangular 
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"plaquettes" is a special feature of the F4 lattice and provides here for the desirable feature 
that (2.5.4) couples only nearest neighbors. 

This feature is desired for Monte Carlo simulations because these do not yield direct 
measurements of physical cutoff effects and one usually estimates the latter by looking at 
the range of interactions in the ultraviolet expressed in terms of the measured estimates 
for the Higgs mass. This range is simply the lattice spacing when the action couples only 
nearest neighbors. It is independent of the values of the couplings in the action. If the 
action contained both nearest neighbor and next nearest neighbor couplings it would be less 
clear what to view as the range and it would be quite reasonable to assume that the range 
does depend this time on the couplings in the action. In the pure nearest neighbor case 
we still don't know, from Monte Carlo alone, how large the observable cutoff effects really 
are, but we can expect with reasonable confidence that variations of the bare couplings 
affect the cutoff effects only through variations in the Higgs mass as measured in lattice 
units. However, when there is a mixture of nearest neighbor and next nearest neighbor 
terms in the action one may suspect that what one calls the "lattice unit" changes when 
the couplings are varied. 

To take the continuum limit we rescale the field = y/6N{PQ + Pi + /?2)^ (we only 
consider the region Pq + Pi+ P2 > The combination ^-^^^qi^qr^p is then analoguous 

to 62 and the combination -^jj^^^^^^^^ analoguous to bi. The first term in (2.5.4) gives 
the usual kinetic term, g{p), on the F4 lattice, 

9(p) = IJ2 [2-cos(p^ + p^)-cos(p^-P^.)] =p'^-^{pY + 0{p^) . (2.5.5) 

The p'^ part corresponds to a negative 60 (2.5.1) which, with the field redefinition (2.5.2), 
generates a positive contribution to bi. Thus the naive nonlinear nearest neighbor F4 lattice 
action corresponds effectively to a continuum action with a positive 61-type term. 

For hypercubic lattices the inclusion of next nearest neighbor couplings is unavoidable 
if one wishes to obtain (2.5.1) in the long wavelength limit. The argument against next 
nearest neighbor terms holds only within the limitations of Monte Carlo methods and is 
therefore relevant to the physical case of N = A. However, at AT = 00 "mixed" actions are 
as useful as "pure" ones because we have at our disposal means to directly evaluate the 
observable cutoff effects. Therefore we shall also consider "mixed" lattice actions. 

As a matter of fact the large N analysis of the model (2.5.4) is quite complicated while 
certain "mixed" models are easier. The reason is that at large N one has to introduce 
auxiliary fields to decouple the terms in the action that are quartic in the fields. One 
would need an auxiliary field for each of the "positive" bonds emanating from a site (12 on 
an F4 lattice) and additional auxiliaries for enforcing the nonlinear constraint. This large 



16 



number of coupled fields makes an analytical treatment cumbersome on the F4 lattice. 
In Pauli-Villars rcgularization one also has quite a few auxiliary fields but most can be 
easily decoupled exploiting Euclidean 0(4) invariance. On a lattice the symmetry is only 
a discrete subgroup of 0(4), and while this might be sufficient, we chose to avoid carrying 
out this exercise. The Pauli-Villars analysis will teach us that all the extra auxiliaries 
are irrelevant to the order in inverse cutoff that we are interested in. But to this order, 
the F4 lattice is essentially 0(4) invariant so the same conclusion should hold on it too. 
Therefore it seems somewhat a waste of effort to struggle to fully solve (to infinite order in 
the inverse cutoff) the pure nearest neighbor action. Moreover, this "exact solution" only 
goes as far as giving us expressions involving some lattice momentum integrals, and, in 
practice, one still has to evaluate the latter by numerical means making it again necessary 
to approximate by ignoring terms that have no effect to leading order in the inverse cutoff. 
To be sure, in Monte Carlo simulations, at A/" = 4, we are advocating the use of the 
"triangle" action (2.5.4). 

We have argued earlier that for Pauli-Villars rcgularization the 62 term plays no role 
to the order in the inverse cutoff considered. We assume that this would also happen on 
the lattice and set /?i to 0. We have also explained that we would like to avoid having to 
deal with the P2 term in (2.5.4). To preserve the freedom of varying the strength of the 
important couplings in the long wavelength limit we need to replace this term by another 
one that has a similar effect at long wavelength but which can be easily "decoupled" 
preferably employing one aiixiliary field only. We choose an action that does this: 

2 



<x,x'> ^ X 



l=<x,x'> 



(2.5.6) 



The parameter 7 has been introduced because the same formula will work for the hypercu- 
bic lattice also; for the F4 lattice 7 = 3, but for the hypercubic lattice 7 = 1. Geometrically, 
the new term corresponds to a coupling between any four fields that live at the ends of two 
bonds that have a site in common. Before, on the F4 lattice, we only had such coupling if 
the two bonds spanned an angle of sixty degrees at the common site. Now larger angles 
are included and fields living on sites separated by more than a single lattice spacing are 
coupled. 

On a hypercubic lattice it is impossible to write down an action like (2.5.4) involving 
only nearest neighbor terms. The second term in (2.5.4), when considered on a hypercubic 
lattice, becomes in the continuum limit proportional to X]^(5/x^ • d^^)"^ and thus breaks 
Lorentz invariance. We therefore have an additional reason to ignore it as we did before 
for the F4 lattice. The action (2.5.6) considered on a hypercubic lattice reduces, after 
rescaling 4> = \/'^N{(3q + (52)^^ almost to the form (2.5.1) with 62 = and ^ (/3o+/32)^ 
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analoguous to bi. The difference is that instead of the bo term in (2.5.1) we now have a 
Lorentz invariance breaking term (p dfi4'- view of this we have no strong reason for 
considering this particular form of the action, except that it includes, with (52 set to 0, the 
case most thoroughly investigated to date for N = A. 

In naive nonlinear models on the hypercubic lattice the Lorentz invariance breaking 
term can be avoided by using a "Symanzik improved" action. For = 4 the exact 
elimination of Lorentz invariance breaking terms at order inverse cutoff square is not 
an attractive proposition because it would necessitate an impracticable amount of "fine 
tuning". At infinite this problem is less severe and we can do quite well with just 
generalizing the (tree level) Symanzik improved action to also contain a hi like term with 
four fields and four derivatives. We are thus led to also consider the action 

Ssi = - 2N/3o J2 (Pi^) ■ ^(^ + - Y^^i^) ■ ^(^ + 2a*)) 



X,IJ, 



20 



(2.5.7) 



As will be seen later on, this choice of the transcription of the four-derivative term exactly 
eliminates Lorentz breaking terms in the full pion propagator at order inverse cutofi^ square. 



3. PHASE DIAGRAM FOR THE PAULI-VILLARS MODELS. 



Similarly to the linear case we introduce auxiliary fields in (2.5.3) to make the depen- 
dence on 6 bilinear. We obtain a new action 



Si = 



+ 



(3.1) 



where uj^^ = (summation over repeated indices is implied) and o^^j/ = coi/jj, and hence 
the S^i, terms above can be replaced by zero. In more condensed notation we have 



^1 



I r ^-.^ N 1 2 1 



(t)K(j) 



[f3p - Ti^lA - -g20Jfii,ujf^u] 



with 



K = K - dfj^Xd/j, + p- d^uj^pdy 
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(3.2) 



(3.3) 



We separate the zero mode of as in the hnear case in (2.1.7) and integrate out tt and H. 
We obtain 



^ 2 



with 



(3.4) 

a^(p-p') = 0, / p' = . (3.5) 
Jx 

3.1. Saddle point equations and dominating saddles. 
The saddle point equations are 
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As promised, the saddle point equations do not depend on g2- Sometimes they may admit 
several solutions. In these cases we need to find the dominating one. Among the possible 
solutions some may have broken symmetry and some unbroken symmetry all at the same 
couplings. In such regions of the phase diagram the order-disorder transition can become 
discontinuous. The analytic continuation of the continuous transition into these domains 
yields metastable critical regimes that have to be eliminated from our search for the Higss 
mass bound. Therefore, a complete analysis of the competing saddles is necessary. We 
shall ignore possible dominating "end point" contributions because we are pretty sure that 
they will not affect our conclusions regarding the accessible critical regime. As will be clear 
later on, the possibility of first order transitions cutting into the critical regime is realized 
and has physical content teaching us something about the dynamics that is both relevant 
and illuminating for the mass bound issue. 

Prom (3.4) it is clear that we wish to find the solution {Xg, Ps) that minimizes the 
following function of A and p: 

i(\ \ /"i P^(l+P^^) + Ap2 + p ^ 1 2 

^(A,p)=/log ^r—r. + vp-(3p+-giX . (3.1.2) 



The function is completely defined when one adds that f ^ > implies p = and p > 
implies = 0. Both p and v"^ are nonnegative. 

To simplify the analysis we introduce some rescaled variables 



u={l + X){167r'^P)^ , t = p(l + A)-'ir . (3.1.3) 
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ra+l 



It is also useful to factor out a positive constant from and consider from now on only 
the minimization of t/^, defined by 

ij{u,t)=un ^d^log + —tu—-tu— + -g*{u-uy (3.1.4) 

where 

gi = 327r2/?V , = (IQtt'^P)^ • (3.1.5) 

It is easy to check that the saddle point equations are reproduced by setting the derivatives 
of ijj with respect to u and t to zero. This of course had to be true. 

We now split the candidate saddles into two classes according to whether the symmetry 
is broken at the saddle or not. In the symmetric phase = and the equation d^l^/dt = 
can be used to define a function t{u) for < u < uq by 

„„ = frT^l^ = f:.3^)^ . (3.1.6) 



t{u) varies between zero and positive infinity when u goes from uq to 0. 

Let Us and ts be the coordinates of a saddle point in the symmetric phase. Starting 
from 

il^{us,ts) = J -^{u,t{u)) + iP{uo,0) (3.1.7) 

i^{us,ts)= / '\G{u)+g\u-u*)]+i;{uo,0) (3.1.8) 

G{u)=u-^ r -—f^-— . (3.1.9) 

^ ' Jo e+^+e+^(«) 

G is defined in the interval (0, uq) and is monotonically decreasing there. 
In the broken phase the equation 



we derive 



where 



'dt 



= (3.1.10) 

t=0 



together with the requirement v'^ > yields the restriction Us > uq. The function x/j at 
the saddle can be written as 

i^ius, 0) = "^u-n ^ + \g\u - u^ . (3.1.11) 
If we extend the range of the function G from u e (0, uq) to the segment u e [uq, oo) with 



2-« r°° EdE 2-« TT 
G{u)^u~ I ^ ^, = u~ — ^^2^ for u>UQ (3.1.12) 



1 + e nsin(f) 
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we can rewrite (3.1.11) as 

ViK,0)= / \G{u)+g*{u-u*)]+^{uQ,Q) . (3.1.13) 

The advantage of these manipulations is that equations (3.1.8) and (3.1.13) have identical 
forms and can be geometrically interpreted. 

G{u) is somewhat complicated but independent of the continuously varying couplings 
- hence it can be computed once n is given. The dependence on the couplings comes in 
through the parameterization of the straight line —g*{u — u*); different couplings corre- 
spond to different intercepts and slopes of the straight line. All candidate saddles are found 
at intersections between the straight line and the "universal" function G. The quantity to 
be minimized is the signed area bounded by G, the straight line, the line u = uq and by 
the particular intersection point under investigation. When the straight line —g*{u — u*) 
intersects G three times and the two areas between the consecutive intersection points are 
of equal magnitude we are at a symmetry breaking transition point of first order. When 
the straight line intersects G only once at the point {uq, G{uq)) the transition is second 
order if the straight line is steeper than the line representing the tricritical case. The 
tricritical case corresponds to a straight line that goes through (uq, G{uq)) and, in addi- 
tion, is tangent to G there. These cases are illustrated in Fig. (3.1) for n = 3. When the 
straight line does not intersect G there is no translationally invariant saddle and we are in 
a frustrated phase. 

In equations the conditions for criticality are as follows: The tricritical point has 
parameters given by 

9t.c. = ^ tan (-) , = 2^«o (3.1.14) 

and the second order line is described by 



3.2. Physical properties of the phase diagram. 

The details of the complete phase diagram are not essential; the schematic structure is 
presented in Fig. (3.2). The "frustrated region" corresponds to regions where translational 
invariance breaks spontaneously and regular saddle points are either not competitive or 
do not exist. 

It is important to understand the source of the tricritical point: As already mentioned, 
the action can also be viewed as an effective chiral Lagrangian with couplings of order one. 
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Figure 3.1 The solid line is the function G{u) for n = 3 and the straight dotted lines 
represent —g*{u — u*) for three typical transition points. All candidate sad- 
dles are found at intersections between the straight line and the function G;p 
the coupling constants in the action only affect g* and u* . 

The single stable particles are the pions. Two pions will attract if they are in a relative 
zero angular momentum and total isospin singlet state. This is easy to understand: If two 
field configurations corresponding to an approximately localized pion are placed one on top 
of the other and the isospin indices match appropriately, the net state will be closer to the 
vacuum than a state where the pions are far apart. The index matching will be right when 
the total isospin is zero. Hence, soft pions attract in the / = 0, J = state, a well known 
fact. One can change the interaction between the pions only at subleading order in their 
momenta, (even the logarithmic part p^logp"^ is fixed by current algebra), and this is 
the main physical effect produced by varying the couplings (of course the value of F^r/A is 
also changed by the variation of the couplings, but this is an "unphysical" effect). If one 
introduces extra attraction, it is possible that two pions of some nonvanishing momentum 
each (in units of the cutoff), can bind to a zero mass state which will be stable and 
condense. We believe that this is what is happening at the tricritical point and we shall 
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Pion Wave Function Renormalization Constant 




Figure 3.2 Generic phase diagram containing the quahtative features common to all 
models investigated. 

present some evidence later on. As far as the Higgs mass bound goes, we want to make the 
resonance replacing the above bound state as heavy as possible. This will be achieved if 
we introduce as much repulsion as possible between the pions, thus delaying the formation 
of the resonance to higher momenta of the "constituent" pions. We shall see that this 
reasoning is born out by explicit calculations. 

To complete our investigation of the region of the phase diagram we are interested in, 
we must also ascertain the local stability of our saddle points (global stability was checked 
in section (3.1) but we ignored until now the question of local stability). To do this part 



23 



of the analysis we need to compute the small fluctuations around the saddle points. Since 
we are interested only in a particular region, our computations will be restricted to that 
part of the phase diagram. We should also remember that some of the field variables 
are auxiliaries. The integration contour for auxiliary fields that are unphysical can be 
deformed in the stable directions, so stability isn't really an issue there. However, we arc 
not allowed to deform the integration contours for the "physical" fields, because if we did 
that we would easily loose the approximate unitarity we have at low energies. Among the 
auxiliaries only p is "unphysical" as its role was just to impose the fixed length constraint; 
however, uj^p and A are simply related to bilinears of field gradients and should stay real 
fields. In short, we must make sure that the pion, Higgs, uj and A propagators come out 
to be positive in Euclidean momentum space. 



3.3. Pion propagator and ghosts. 



From now on we shall always assume that the couplings are chosen in such a way that 
we are somewhere in the broken phase close to the second order transition. 

We again separate out the zero mode from the field (j) and parametrize the remainder 
by H and tt (see (2.1.7)). Expanding around the saddle we introduce the shifted fields 

5\ = \ — \s , SH = H , Stt = TT , 5uj^y = uj^y , Sp = p . (3.3.1) 

From (3.3) we read off the pion propagator in Fourier space 

xab xab 

< Stt'^Stt^ >= — ^j- = — r^TT-TT ■ (3.3.2) 

K{p^) + XsP^ (1 + As)p2 + (p2)n+l ^ ^ 

There are ghosts in this propagator (poles with residues that are not real positive num- 
bers) . Their presence refiects the fact that Lorentz invariant Pauli-Villars regularization is 
achieved at the expense of exact unitarity in Minkowski space. The ghost poles are located 
on a circle in the complex plane 



\plhost\ = + ^s)^ ■ (3-3.3) 

The condition A,, > —1 is always satisfied. When the energies in a process approach 
(1 + As)2^rA = As one expects violent cutoff effects to set in. Thus, the "physical" cutoff 
scale is coupling dependent and different from the "bare" cutoff A. In retrospect we see 
that it makes more sense to measure our dimensionful quantities not in terms of A as we 
did until now but rather in terms of Ag. When we do the lattice analysis we shall see that 
violations of Euclidean rotational invariance occur in the pion propagator at a distance of 
the order of the lattice spacing and that there is no dependence on the couplings. In our 
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case the correct analogue of the inverse lattice spacing is Ag, not A, and rescaling by As 
here is analogous to the standard practice of setting the lattice spacing equal to unity in 
lattice work. 

We rescale all our fields (including the auxiliaries) and couplings, but keep the old 
notation, 



a^, p d^, p X (1 + As) 



j_ 

2n 



dp, X ^ dp, X X (1 + As)~^ 
H, 7T ^ H, jf X (1 + As)^ 

A ^ A , , 

(3.3.4) 

p p x(l + As)" 

U ^ UJ 

V ^ f X (1 + As)^ 

2 

91,2 91,2 X (1 + As)" . 

The original constraint ^ = now becomes ^ = {l+Xs)~^N/3 and the pion propagator 
becomes 

< S^-Stt' >= ^ ^^^^2 + (^2)n+l] = (1 + Xs)K{p^) ■ ^^-^-^^ 
With standard conventions we define the pion wave function renormalization constant Z^^ 
by 

= ^ . (3.3.6) 

This constant is fixed by the couplings through the saddle point equations. From our 
analysis of the phase diagram we deduce that the range in which is allowed to vary is 

0<Z,<'^-^^ . (3.3.7) 

depends only on g\ and it is useful to invert the relationship viewing Z-,^ as the external 
control parameter restricted only by (3.3.7) 

1 <^ ' ' (3.3.8) 



Ik K{k) 

(3.3.8) is a rewriting of the second saddle point equation and fixes As. The other saddle 
point equation fixes v. We know that ultimately we shall be more interested in the quantity 
jZ-j^ because it is equal to the pion decay constant. So we choose to write the first saddle 
point equation in the form 

~ '^^'^ Jk K{k'^) 
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3.4. Small fluctuations in the broken phase. 



After integrating out the pions from the action (3.2) and taking into account the 
rescalings (3.3.4), we expand around the saddle point. To quadratic order in the small 
fluctuations we obtain 



(2) 



N-1 2 
2 



^ [ SHKSH + VNv [ SpSH + ^ [ {9idX^ + g2Su;l^)+ 

'^TT Jk Jk o Jk 



UpBPPSp + SpBPhX + SpB'iPSuj^y 



+ 



(3.4.1) 



where 



B^\p^) 
BP\p'^) 



Ik 



kK{\p + k)K{lp-k) 
y - 

kK{\p + k)K{\p-k) 
\PliPv - k^kv - \8^v{\p^ - k'^) 



Ik 



k K{\p^k)K{\p-k) 
_ j [\pi,Pv - k^ky - |V(|p^ - fc^)](y - k'^) 
Jk K{\p^k)K{\p-k) 

[\pnPu - k^ky - \5^j,y{\p^ - k'^)][\p^,'Pu' - k^tkyi - y^'y'iy - k"^)] 



K{\p + k)K{lp-k) 

(3.4.2) 

To compute the various propagators we would have to invert a 12 x 12 matrix. To 
simplify this task we exploit rotational invariance to block diagonalize the matrix. We 
decompose the field Su^u in a spin-zero, spin-one and spin-two field. To do this we 
introduce Euclidean polarization vectors W^(p), where j = 1,2, 3, 

Wl{p)W^{p) = 5^^ , p^wlip) = . (3.4.3) 

We now define new u fields 

SuJ^iy{p) 

I 11 

Wi{p)wi{p)5ojjk{p) + T^[p^,Wl{p)+pyWl{p)]5ojj{p) + -^{pf^pu - -Sf,r,p^]Su;{p) . 

Z\p\ p 4 

(3.4.4) 
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The spin-two field Sujj}. is symmetric and traceless. Because of Lorentz invariance the 
12 X 12 matrix must now spht up into three blocks: a 4 x 4 block involving the scalar fields 
5H, 5X, 5p and Su, a diagonal 3x3 block for the vector field Sujj and a diagonal 5x5 
block for the tensor field Sajjj^. This decomposition is exact in Pauli-Villars regularization 
because of the preservation of rotational invariance; on a lattice this decoupling would not 
be exact, but, as long as one is ultimately going to expand in the inverse cutoff only to 
leading and subleading order, a similar simplification should occur for F4 lattice actions 
that don't break the lattice symmetries. In terms of the new 5oj fields the quadratic action 
becomes, ignoring order 1/N corrections. 



= Ar I 8HK5H + yfNv [ dpSH + —Z^ 
2^7r Jk Jk 2 

/ SXil-B^^ + -T^gi)SX + SXE^'^Su + \5pBPP5p + SpRPHx + SpB^'PSco + ^Z^ 
Jk I 2 4Z^ 2 J 2 

f I 1 3 11 11' 

/ Su{-B'^'^ + Y^g2)Suj + 6u;j{-B'^'^ + ^g2)5ujj + 5ujjk{-Bf^ + -^g2)8ojjk 

(3.4.5) 

The B "bubbles" in (3.4.5) are projections of the B bubbles in (3.4.2) via (3.4.4). Their 
explicit forms can be read off equation (3.4.8) below. Introducing the four component fields 
il)js^^ representing Sp, 5X and Suj for A = 1,2,3,4, respectively, we rewrite eq. (3.4.5) 
as 



4'^ = 1 f ^aMab^B + 1 f [MvJ^i^oojf + Mt J^i^oojkf] . (3.4.6) 
Jk Jk ■ 



The variables M are functions of the couplings and momentum square; they can all be 
expressed in terms of six elementary "bubble" integrals In,m{p^)^ m + n <2: 
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^AB is symmetric with non-vanishing entries 
Ml2 = "//Vv 

M23 = InzI (^^/o,o - /i,o) 

, 2 

-2\2 



M44 = -WZ, - ^55-^0,0 + ^/o,i - h.2 - -^hfi + j/l.l - ^hfi j 



3.5. Spontaneous breakdown of 
space-time invariances and local stability. 

We are finally in position to check local stability for the vector and tensor fields. The 
corresponding matrix entries Mj- and My first vanish at zero momentum when g2 decreases 
to a critical value g2c given by 



At g2 = g2c ^Iso the entry M44 vanishes at zero momentum. In all our subsequent cal- 
culations we shall assume ^2 to be safely larger than this critical value so that even for 
negative (but small in absolute magnitude relative to unity) values of these entries are 
of order unity and positive. Physically, we wish to keep the masses of the vector and the 
tensor of the order of the cutoff. Note that both fields are isoscalars. 

On a lattice a related phenomenon would be that a non-translational invariant saddle 
takes over. In Pauli-Villars regularization when translations get broken so do rotations 
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and we obtain both a massless vector and a massless tensor. On a lattice rotations aren't 
a continuous symmetry and we would not expect a "tensor" like particle to also become 
massless. 

Thus, although g2 doesn't appear in the phase diagram there is some bound on it too. 
Once this bound is satisfied we have no evidence for any other source of local instability 
and we believe that the 4x4 matrix M will be positive definite for all Euclidean momenta. 
We have not tested this fully, but from the behavior at low momenta, which we did test 
(see below), it seems very safe to assume that no local instability is hiding at high momenta 
and that all critical regions that we shall henceforth be interested in are indeed accessible. 



3.6. Spontaneous breakdown of scale invariance at the tricritical point. 



The limitation on in (3.3.7) is turned by eq. (3.3.8) into a limitation on gi 

1 27rn2(n-2)sin^ 
-oo < — < -2 ^ . 3.6.1 

The tricritical point is at gi = gi^t.c. where 



9l,t.c. = - — ^ — ll . 27: ■ (3.6.2) 
27rn^(n-2)sm^ 

Note that the rescalings (3.3.4) have changed the definition of gi so that (3.6.2) differs 

2 

from the combined effect of (3.1.5), (3.1.6) and (3.1.14) by a factor of . 

When gi approaches gi^t.c. it is easy to check that the matrix element M33 from (3.4.8) 
vanishes at — 0. This implies that det M{p'^ = 0) = 0, because, at — also Mn and 
M34 vanish and then the first and third row of M are proportional to each other. The zero 
eigenvector is a particular linear combination of 6H and 6X. The field A is related by the 
equations of motion to (5^0)^; this is obvious from eq. (3.1). A pole at zero momentum 
in 5X is therefore likely to be interpretable as a dilaton. Therefore, at gi = gi^c the Higgs 
particle becomes massless and also plays the role of a dilaton. This behavior is similar to 
the one discovered by Bardeen, Moshe and Bander in three dimensions with a bare action 
that was renormalizable [17]. 



4. HIGGS MASS BOUND AND CUTOFF EFFECTS 
WITH PAULI-VILLARS REGULARIZATION. 



We are now ready to address the main problem: How large can one make the ratio 
m^//7r while keeping cutoff effects small? We shall carry out our analysis in two stages. 
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At the first stage an approximate but simple calculation will tell us what to expect. At 
the second stage we shall perform a complete analysis. 

First we need to make an observation that will simplify matters for both stages: The 
Higgs mass will be obtained from the matrix of (3.4.8). We are only interested in 

small values of /tt (i.e., F-^- in units of As) and for such values the Higgs resonance will 
appear at small complex value of (also in units of Ag) because the physical coupling, 
even if large, is a finite number. Therefore, it is consistent to expand the entries Mj^j^ijP") 
in powers and logarithms of p^ . In ordinary, rcnormalized contimmm field theory (i.e. 
when calculating only universal quantities), we would stop the expansion at leading order. 
Here we shall be going to one order higher. The observation we wish to make is that 
to first subleading order the scalar 5ul> decouples and the Higgs mass and width can be 
obtained from the 3x3 upper left corner of the 4x4 scalar block of M. We shall denote 
this submatrix by M''. 

In addition, 5u3 contributes to tt tt scattering, given by the exchange of the scalar 
auxiliary fields, only at orders that we shall be neglecting. Almost all the information 
needed for calculating tt tt scattering to first subleading order in the inverse cutoff is 
contained in M^. Fluctuations in the approximately decoupled low momentum components 
of 5u} are under control because, as we have seen, local stability requirements for the 
vector and tensor low momentum excitations stabilize also the scalar 5uj. Thus, finally, 
the dependence on the precise value of g2 also disappears from all the processes we shall 
be interested in. 

Once we decided to go only to first subleading order we have to focus on M'^ and from 
(3.4.7) and (3.4.8) we see that we only need 3 out of the 6 "bubble" integrals in (3.4.7), 
namely q fo^' n <2. For arbitrary n it would be very inconvenient to try to work out 
closed forms for the complete integrals as we did in the linear case. Moreover, unlike in 
the linear case, we have already dropped some higher order terms when we reduced our 
problem to so we should be consistent and keep only the leading and subleading terms 
in the external momentum for the three bubbles that we need. Some of the technical 
details of the relevant computations are sketched in Appendix A. The results are 

/0,0(p') = T7^[-logp2 + 1 - 1 - (1 - \)-^p''] + ■■■ 
16-71^ n 12 sm^ 

hfi{P^) = TTT^ 1 - -)^^ + Tr log/ + / + • • • 4.1 

' 167r^ n nsm^ 4 12n ^ ' 

2^ 1 r (»-2)7r (n-l)(n + 4)7r 2. , 
IGtt^ n^sm^ 12n^sm| 

Note that the limit n — > oo can be taken on the leading terms but not on the subleading 
ones. This refiects the fact that with a sharp momentum cutoff the first correction is 
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suppressed by \p\ and not by p which means that non-local operators come into play. 
Also, there is a slight error in the equation for the simplest bubble, /q^q? [S^IO]: at 
infinite n we get to leading order in p^, /o,o(p^) — \^[~^^^p'^ + ^i' while in [9,10] 
Iq,q{p^) = iQ^2 [~ logp^ + 2] is used. This has no significant effect on the physical numbers 
obtained in these papers. As a matter of fact, many workers adopt a convention where 
the "physical" cutoff Kj^ (the "Landau pole" ) is defined as the point where the asymptotic 
expansion of /q^o, las — i^f2[~logP^ + c] obviously breaks down, namely, /as(A|^) = 0. 
With such a convention, the particular value of the constant c does not need to be known, 
and all dependence on the cutoff scheme disappears. This convention is of course quite 
arbitrary but not totally unreasonable; it does explain to a certain extent the effective 
"universality" of the triviality bound. 



4.1. Approximate calculation. 



It is much easier to do calculations for weak couplings. Ultimately we wish to find out 
how strong the physical coupling can become without distorting the theory too much. This 
is achieved by limiting the cutoff effects. If we make this limitation extremely stringent the 
cutoff is very high and because of triviality we are forced to weak physical couplings where 
it is easy to calculate. Even for very stringent limitations on the cutoff effects there will be 
a dependence on the bare couplings and in some region of the space of bare couplings the 
physical coupling will be allowed to be larger (while still very small in absolute magnitude) 
than in other regions. It is very reasonable to assume that there is some smoothness in 
the dependence on the higher dimensional operators that is induced by the variation in 
the bare couplings and, therefore, the region we shall be interested in is the one where we 
expect the largest possible coupling even when the cutoff effects are less stringently limited. 
In particular, we already know that in practice the physical coupling cannot be made large 
in a sense that would invalidate perturbation theory completely. The major result to 
date is that even for relatively small physical couplings the cutoff effects become sizable 
disallowing further increases. Therefore, our decision to first work where the physical 
coupling is small is not expected to lead us astray. 

When the physical coupling is small the width of the Higgs particle is small too and, 
to simplify the formulae, we shall ignore width effects and concentrate on the quantity mji 
defined by 

Re[detM^(p2)]^,_^^ = . (4.1.1) 

mji is related to, but not equal to, m^. It can become quite different from mu even for 
moderate physical couplings but the relationship between the two is monotonic and, for 
the purpose of identifying the right region in the action space, mji is perfectly adequate. 
Similar considerations were presented in our study of the linear models. 
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In our approximate analysis all the subleading effects are small but still hard to cal- 
culate. Since they are small, they do not affect the numerical value of the leading order 
answer.* We could do without calculating them explicitly if we had an independent means 
to tell when a given case will have larger cutoff effects than another case. The easiest is to 
use the pion propagator, because it has a simple explicit form to all orders in the inverse 
cutoff (see eq. (3.3.5)). In units of Ag we know that ghosts will appear for energies of order 
one and this is independent of the bare couplings. Therefore, making niji in units of Ag as 
small as possible while keeping the ratio mji/ fj^ constant will diminish the cutoff eflFects 
and identify the right region in the space of bare actions where a larger triviality bound 
ought to be found. To be sure, our approximation isn't quite consistent logically because 
the cutoff effects on the pion propagator are suppressed by more than one power of and 
hence do not reflect the dimension six operators. We still think that our approximation is 
useful because, if one does not "fine tune" in the sense explained before, once the cutoff 
effects induced by the dimension 6 operators become sizable (of the order of 10 percent) 
all the higher dimensional operators also quickly become important. Therefore, the pion 
propagator does sense the point where cutoff effects turn on in an overwhelming way. Our 
argument is also helped by the fact that our more complete analysis confirms the findings 
we shall present here. We should point out that most lattice work was essentially equiva- 
lent to viewing the inverse lattice spacing as ~ Ag/yr. This point of view is not logically 
superior to the one we adopt in this section but nevertheless produces quite reasonable 
results when compared to the outcome of a more sophisticated analysis. 

By general arguments we know that 



m| ^ C2(Z^) exp[-967rV^ii] (4.1.2) 



where 



2 



^i? = , = ^ , (4.1.3) 

and we have already used the fact that any desired value of f-^ can be obtained for any 
by tuning Therefore the constant C depends only on Z^j- (which has been traded for 
gi in (3.3.8)). All that is left to compute is this constant and for this we only need the 
leading terms in eq. (4.1). We plug in the needed results from (4.1) in (3.4.8) and then in 
(4.1.1) ending up with 



C{Zt^) = exp 



n — 1 / TT A/ -7 \ 

^ [ + ntan(7r/n)^^ 



(4.1.4) 



l-Z^(n-2)/(2(n-l)) 



* 



A quantitative feeling may be obtained from the analogous situation in the linear case: see 
the expansion in eq. (2.4.2). 
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Since, as explained above, we want to find the range of Z-,^ where, for a fixed ratio 
mjil f-!^, mji is minimal, we need to minimize C. We must remember that Z^^- is restricted 
by (3.3.7). The fact that Z-^ has to be positive is no surprise and we now see that, at the 
other end of the interval, where Zjy = 2{n — l)/(n — 2), C diverges showing that the critical 
region shrinks to zero there. This is where the tricritical point is located and our analysis 
confirms our physical arguments that indicated that the bound would be largest when one 
stays as far as possible from the tricritical point. 

As far as the numerical validity of the approximation is concerned we can get some 
feeling by looking at the higher order (neglected here) terms in (4.1) or at the expansion 
in (2.4.2) for an analogous problem. An accurate evaluation will be given later on in 
subsection (4.6). The numerical values one gets from (4.1.2) are, for example, good to one 
percent if mji corresponds to a physical Higgs mass no larger than 0.800 TeV, Z-j^ <1 and 
n < 100. For larger values of n the accuracy deteriorates somewhat. 



4.2. Some numbers. 



To get a feeling for the numbers involved let us take = 4 and ask what change in C 
will induce a noticeable change in away from the value = ^rys ~ 2.5 corresponding 
more or less to the present bound. Maintaining the cutoff effects at a fixed magnitude with 
5{mji) = we get 

-aogC = = ^^(^) - 6 i&) = 2mi^ (in TeV] . (4.2.1) 



To increase by 0.100 TeV we need to decrease C by a factor of e ~ 7.4 Between 
= and Z-j^ = 1 little variation is induced, but when approaches 2^^^, will 
decrease sharply. In summary, for almost all large enough n's the region < Z^^ < 1 
seems to give reasonable estimates for the bound. In terms of gi this region corresponds to 
gi < 0, showing that the parameterization of the space of actions we chose is a reasonable 
one and lending support to our physical understanding of the mechanism affecting the 
bound. 

When = 1 = by (3.3.8) and the model is the simplest nonlinear model possible. 
While this is not the most optimal place to look for the bound, it is reasonably close 
to it. Therefore, even the simplest model, and even setting n = oo, as we discussed 
above, will give us reasonable and quite high estimates for the bound. On the lattice the 
situation is not exactly the same, however, and this is the main reason for present estimates 
being quite significantly lower. A lattice action of the naive type has generically, in the 
notation of eq. (2.5.1), a nonvanishing negative parameter bQ in the derivative expansion. 
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To approximately translate this into a Pauli-Villars regularization one must redefine the 
fields according to (2.5.2) and this induces an effective gi coupling (and some g2 coupling 
also, but this has no effect) with a positive sign. Hence naive lattice actions are somewhat 
outside the region where reasonable estimates for the bound can be obtained and, because 
of this, the lattice bound is expected to increase. We shall see that the increase is not 
negligible, but not very dramatic either. 

4.3. Computing leading cutoff effects in 1/N. 

We have reached the point that we wish to do a complete calculation. In this subsection 
we explain in detail the calculation conceptually. 

The physical scale is set by the unitless /t^. Therefore, we consider /^r as an exact 
function of the bare parameters. By definition there are no "cutoff effects" in fi^. Similarly 
we define another "exact" quantity, the coupling g = 3m|^//^. mjj is the exact real part 
of the pole of the amplitude for / = 0, J = w- tt scattering when continued to the second 
sheet below the physical cut. mjj is a function of the bare couplings with no "cutoff 
effects" by definition. 

Consider now some new physical quantity. We make it dimensionless by extracting 
the appropriate power of f-^ and denote it by P. P may depend on momenta Qi which we 
measure also in units of /t^: qi = rifj^. Let the bare action depend on n bare parameters. 
We imagine changing variables to /tt, g and n — 2 remaining parameters p and consider P 
as a function of them: P = P{ri, r2, . . . ; /tt, g',p)- The change of variables holds in some 
neighborhood of a particular point in the broken phase that we wish to investigate. We 
keep p fixed and carry out a double expansion of P in and g without paying attention 
to whether the change of variables is one to one also in the region where this expansion is 
sensible. Note that the dependence on /tt is both explicit and implicit via the momenta. 

If we expand P in g at fixed r^, f-^- and p renormalizability and universality tell us that 
to any finite order in g the limit /tt — exists and is independent of p. The summation of 
all orders in g of the /tt = terms is ambiguous in the full theory because the series is badly 
behaved. However, in our case we expand P in 1/N after the appropriate rescaling and 
redefinition of ^. To leading order in 1/N P has a nontrivial limit as /t^ — > even without 
expanding in ^. It is unlikely that this continues to be the case at subleading orders in 
1/N because one cannot replace the leading order (in 1/N) propagators in higher order (in 
1/N) terms by their leading asymptotic expressions in f-^ as this will induce divergences 
in the momentum integration at the positions of the "Landau poles" . But at leading order 
we are lucky and this simplifies our task considerably. 

The above may look like a contradiction to triviality, because we just argued that at 
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leading order in 1/N we can obtain a nontrivial continuum limit. The catch is that the 
change of variables from the bare parameters to f^^, g and p excludes precisely the region 
where the continuum limit is taken. There are no values of the bare parameters that are 
physically acceptable and can maintain a nonvanishing coupling g while /tt is taken to 
zero.* However, we can carry out our change of variables away from the critical line in 
the broken phase and work out the expansion by analytically continuing (even without 
noticing) into the unaccessible region.^ In this way we can, at least at leading order in 
1/N, disentangle the issue of triviality from the issue of summability of the renormalized 
perturbation theory. 

Once we have the fw^O limit, PN=ooiTlif2i • • • ! 0; S'); we can expand in the inverse 
cutoff 

-PAr=oo(n, ?'2, • • • ; fl, g; p) - PN=ooirh '^2, • • • ; 0, ^) = 



APAr=oo(n, r2, . . . ; fl g; p) + 0{fi log^U) 

where T is some finite number. We define the leading cutoff effect by 



(4.3.1) 



PN=oo{n,r2,---;0,g) 



Sp{ri,r2,...;U,g;p) ^ — j:^—: rTTTN — • (4.3.2) 



In practice we always work on a one dimensional line connecting the given point in the 
broken phase where we wish to calculate the cutoff effects to some point on the critical 
line. We change variables in our parameter space (and this change is globally well defined 
in a large region that includes the critical manifold) from the original n parameters to 
n — 1 parameters p and one parameter ^. ^ can be traded in a one to one fashion for g. 
Varying p and g we span the whole region we are interested in. We organize our work by 
fixing p and varying g along the line. Our purpose is to find the region in p where g can 
be maximized under the restriction that the cutoff effects on some physical observable(s) 
do not exceed a given amount. Along the line /tt is monotonically dependent on g and the 
cutoff effect is measured by 5 where 

Sp{ri,r2,...;g-p) =Sp{ri,r2,...-J^{g),g]p) . (4.3.3) 



In practice, we cannot compute the function g exactly in terms of the bare parameters 
and cannot carry out the change of variables explicitly. But, we can compute S exactly, 

* This can be seen as follows: For small g and fixed Z^^, mjj satisfies an equation that is 
similar to (4.1.2) and for all acceptable Z-,^ one has C{Zt^) > Cmin > 0- Replacing mjj on 
the left hand side of the equation by /t^ and g we see that if g is fixed /tt is bounded from 
below away from zero. 

t There are claims that if this is really done with care one will discover ambiguities and these 
are a direct refiection of the singularity structure of P in the Borel variable conjugate to 

g [18]. 
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because errors of order log ^/^^ or higher in g have no effect and it is practical to compute 
fn{g) neglecting order log^f^^ contributions. In plots we shall always show ^ as a function 
of Mjj = \/g/3 X 0.123 TeV in order to be able to immediately read off the Higgs mass 
bound in TeV. 

There is a fundamental difference between our way of calculating cutoff effects here 
and the one employed previously in the triviality context. Here, the single approximation 
in computing S is the expansion in 1/N. In other computations, at = 4, one uses the 
loop expansion. In terms of the bare coupling, the loop expansion is a finite rearrangement 
of the perturbative expansion. An answer correct to / loops will be also correct to order 
/ in the bare coupling. This is still true when the series is reexpressed in terms of an 
appropriately defined renormalized coupling, g. The series for S in g has the following 
general structure: 

oo n 

S = g^fl E ^" E^n,/(log/^) (4.3.4) 

n=0 ;=o 

where Vji j is a polynomial of degree /, / is the number of loops, and T is some number. 
Along a line in the space of bare actions that touches the critical surface at its end we 
have /tt ~ exp[— 487r^/(7] (at = oo) and one needs to go to an infinite number of loops to 
get 5 to a reasonable accuracy, even when g is small. Therefore, estimates of cutoff effects 
obtained by truncating the loop expansion are not under good control. In Appendix B we 
give an explicit example of this problem. 

4.4. Leading corrections to the width to mass ratio. 

The Higgs resonance is a complex root of det M''(p^) = on the second sheet 

detM'^(-(mjy-^7if)^) = • (4.4.1) 

What is meant by "second sheet" is that the branch of the logarithm is chosen so that, for 
< Ih/^H « 1: 

\og{-{7nH - '-IH?) ^ log(r4) -i^-iTT . (4.4.2) 

In practice it will turn out that we have to restrict our attention to the region where the 
width does not exceed the mass by much. 

With rescaled bubble integrals 



Bjip'^) = 167r2/^-,o(p2) 
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(4.4.3) 



the explicit form of eq. (4.4.1), for = —{mjj — ^'Jh)'^ is 
2?^M = _B„(„2) ||p^Bo(p^)-Bi(p^)]^ 

Our problem is to solve this complex equation for the real unknowns rajj and as a 
function of to first subleading order (logarithms are counted as order one) for small 
f^. We no longer assume weak coupling, so the ratio mjj/ f-,^ is taken to be of order one 
too (up to logarithms). 

We are working along lines of constant Z-j^ which plays the role of the set of parameters 
p in section (4.3). It is useful to introduce a convenient parameterization of the line which 
is chosen so that the equation simplifies. On the line we have three real unknowns, all small 
relative to unity: /tt, rriu and 7/^. (4.4.4) gives us two relations among the unknowns and 
the remaining free parameter describes the line. We choose the free parameter to be an 
angle 6 defined by 

irriH - ^ifH? = /^^ exp[-i6'] (4.4.5) 

where /x^ is real and positive. We shall be interested only in that portion of the line where 
the mapping between 9 and the coupling 

2 
Jit 

is one to one (in practice this limits the magnitude of g but cutofi^ efifects become large 
before this limit is reached). This portion of the line includes the critical point at one of 
its ends. 

The two real equations (4.4.4) have to be solved now for fi^ and g in terms of 9. We 
shall obtain the solution by first solving to leading order and then perturbing around it 
with the subleading terms. The leading order solution is 

9o + 7r 



g = 3^ (4.4.6) 



fiO = C{Ztj:) exp 



2tan(^o) 



967r2cos2(%)sin(^o) 
90= — 



(4.4.7) 



(^O + tt) 

It is easy to check that, as ^0 we get for mu the same equation we wrote down for 
mji in (4.1.2). Note that ^o(^o) is bounded in the interval of interest, ^ <9q<'k. 



The physical quantity, P(;0, (7) = we are interested in is dimensionless, and has 
no external momenta (we suppress the subscript N = 00 that we used in section (4.3)). 
To leading order in the cutoff we have 

P(;0,^0) = 2tan(^) (4.4.8) 
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where is given in terms of ^ro by the inverse of the second equation in (4.4.7). 

To calculate to next order we treat all three variables 6, g and //^ on equal footing. 
These variables are defined by (4.4.5) and (4.4.6). At leading order we had 6 = 9q, fi = fiQ 
and g = go given by (4.4.7) with $q as a free parameter, a coordinate along the constant 

line. If we treat all three parameters on an equal footing we ought to also allow for the 
subleading terms to change the coordinate along the line. We set 

/^^ = ^0 + 

6^00 + 50 (4.4.9) 

g = go + Sg ■ 

The equations (4.4.4), when expanded to subleading order give us two relations among 
the three 5's above. The needed third relation comes from noting that the quantity 
P{'i fn^ 9'i Zn) in (4.3.1) is evaluated at the same g as the quantity P(; 0, g) there. Therefore 
we have 

6g = . (4.4.10) 

Carrying out the algebra we obtain 

S Th iao; Zn) = r T^T T (4.4.11) 



M 



H 



1 + (tt + ^o) [tan (^%) - cot (6*0) 



where, on the right hand side, $q and /xq are functions of go defined in (4.4.7). gQ is a 
free parameter and the subscript has no particular meaning any more. The function 
D{Zt^, n) is given by 

nf7 ^ 7r(n2-l) ^ (n + 3)cos^ ^^^ ^(^ _ _ 2) cos^ ^ ^ 

^(^-"^ = 12;?^ + — 6 — + — i^;;^- — ^ • 

(4.4.12) 

The function ({ZTj-,n) was defined in (4.1.4), only now we have written out explicitly its 
dependence on n. 

The overall structure of eq. (4.4.11) is very simple. The correction has factorized into a 
function that carries all the dependence on the coupling g and is universal times a function 
that carries all the information about the cutoff, through the dependence on Z^j- and n, 
namely C'^{Zt^, n)D{ZT^, n). It is amusing to note that the dependence on Z-,^ always comes 
in through the combination that appears in the definition of 

This factorization is more powerful than what one would have expected on the basis 
of the Symanzik improvement formulae, as they are usually presented, in several respects. 
First, we seem to have only a single operator insertion because if we had a linear superpo- 
sition of several operator insertions nothing special should be evident when one looks at a 
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single momentum independent observable as we are. Second, the coefficient of the operator 
is independent of the renormalized coupling. Third, the logarithms in the subleading terms 
have been summed in closed form. To be sure, we are looking at a physical operator, not 
an arbitrary correlation function, so the operator counting should be different. Also, we 
are working at = cxo where additional simplifications should occur. Our result could 
be explained if it were true that, at = oo, all cutoff effects in on-shell dimensionless 
physical quantities (that are functions of dimensionless momenta) are given by an effective 
renormalized action, 

Seff^SR + ce^^[-mT:'^/g]0 , (4.4.13) 

where O is a renormalized operator, c is a ^ independent free parameter containing all 
the non-universal information and Sji is describing the usual universal part of physical 
observables with the unit of energy set by /tt = 1. In terms of general RG reasoning this 
representation of the S^j j is not unreasonable if one accepts that at A^ = cxo the number 
of available independent operators decreases. We shall see that this would also explain our 
results for the tt-tt scattering amplitude. 

D{ZT^,n) is a quadratic polynomial in the variable ( and has two negative roots. It 
turns out that one of the roots can be realized in the allowed range of Z-,^. For n = 3 the 
root is at Z^^ 0.52 while when n — > cxo the root is realized by Z^^ fa We see that 
for all n the roots are for Z^^ between 1 and 0. If we set Z-^- to be at the appropriate root 
for the value of n chosen we have achieved "improvement" in the sense of Symanzik in 
that the leading cutoff correction has been made to vanish. This would be the type of fine 
tuning that we excluded because, as far as we know, it is unreasonable to expect the more 
complete theory in which the minimal standard model is embedded (assuming that this 
is the situation in nature) to conspire to achieve such a more refined decoupling. We are 
entitled however to use a reasonable range, say < Zj^ < 0.5 and n < 100 to estimate the 
bound. 

We find that for Mjj < 0.820 TeV the cutoff effect on the width to mass ratio is less 
than half a percent, as will be shown in Fig. 7.3d. At such a high mass at infinite A^ the 
system is strongly interacting. 

4.5. Leading corrections to tt-tt scattering 

We now repeat the analysis of the previous section for tt-tt scattering. We are con- 
sidering the process tt'^(I) + tt''(2) — > tt'^(3) + tt^(4) where a, b, c, d are isospin indices and 
(1-4) denote momenta. We work in the center of mass frame and wish to compute the 
invariant amplitude A{s,t,u) defined in (2.3.9) to first non-vanishing order in 
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To leading order 'm.l/N and up to and including subleading order in the inverse cutoff, 
t, u) is dominated by the exchange of scalar-isoscalar fields and therefore is just a 
function of s, A{s). The differential cross-section with identical isospin indices for the 
in-coming pions, in the center of mass frame, is isotropic and given by 

/ da\ -/V" , , , , , 9 

37^ =w:r^Ms) 4.5.1 

where s is the square of the center of mass energy, henceforth denoted by W, and we are 
working to leading order in 1/N. 

In fact, keeping only the leading order in alone, at any finite cutoff, is enough 
to exclude the contribution of exchanges of the spin-one field, Suj, by Bose symmetry. 
However, the exchange of the spin-two field, dujj^, does come in, only at one order higher 
in the inverse cutoff than the order we are calculating to. In the case of a linear action 
we had a similar simplification. An immediate consequence of this is that we shall see no 
higher spin resonance in our computations (a p-like particle for example). To see an even 
spin resonance it might be sufficient to just keep all orders in the cutoff, but to see an odd 
spin resonance one must go to a higher order in 1/N. 

To evaluate A{s) and see where the observations in the previous paragraph come from, 
we need to write down the pion interactions. We go back to equation (3.2) and introduce 
the pion fields and the Higgs field by eq. (2.1.7): 

If N f r N f 1 1 

Si = - [tt^tt + HKH] - p -^v pH-- [(3p- -giX^ - -g20J^^uu;^u] • 

*/ X *J X J X J X 

(4.5.2) 

Here K is defined in (3.3) as K = K — d^Xd^ + p — d^w^i/di,- We expand around the saddle 
point, carry out the rescalings (3.3.4), but do not integrate out SH and Stt yet. There 
are three types of vertices involving pions, all involve two pion fields with isospin indices 
contracted and a third field, which is 6p, 6X and 5a;^j/, respectively. At leading order in 
1/A^, A{s) is dominated by diagrams containing two vertices of the above type and the full 
propagators (to leading order in 1/A^) of the ifj fields ((3.4.6)) and of Sujjjf ((3.4.4)). 

The expansion in the inverse cutoff is an expansion in where, up to logarithms, 
rrijj and all momenta Pj are of order unity in We are still assuming that g2 is safely 
larger than g2c ((3-5.1)) and therefore the spin-two field has a "mass" of order imity. We 
also know that ultimately we shall continue analytically to Minkowski space with on-shell 
external pions, p| = 0. Under these conditions we first carry out an order of magnitude 
estimation of the vertices and propagators that can be read off (4.5.2). For this purpose 
it is better not to think in terms of the full propagator, but use multiple insertions of 
the two point vertices in (4.5.2), which amounts to the same thing. It is then easy to see 
that the only kind of diagrams that contribute to leading order have the external pions 
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connected to two TTirSp vertices. The internal line has multiple insertions of 5pSH vertices 
and of the various tttt bubbles. When we decide to also keep the first subleading terms we 
have to add diagrams where two of the external pions are still coupled by a rnrp vertex 
but the other two external pions couple to a ttttA vertex. The internal lines in these new 
diagrams have essentially the same structure as before. The conclusion is that, in terms 
of full propagators, we have, with Euclidean q — Pi + P2, but keeping only terms that will 
contribute to the on-shell amplitude, 



The full propagator < dpSp > (g ) is needed to leading and subleading order, while the full 
propagator < dpdX > {cp') is needed only to leading order. These propagators are obtained 
from the inverse of the matrix M (3.4.6,3.4.7,3.4.8). Some more order of magnitude 
estimates show that the 4x4 matrix M can be replaced, to the order we are interested 
in, by the 3x3 matrix M'' defined at the beginning of section 4. 

In the notation of section (4.3), the dimensionless quantity we are interested in here is 
chosen to be 



where r = W/F-,^. One expects a bump in P when W = M^j and the width is small. It 
therefore makes more sense to replace r by another dimensionless number, R, so that in R 
the bump would be independent of g when the width of the Higgs is small. We therefore 
define 



effect we are required to keep r and g fixed and this is equivalent to keeping R and g fixed. 
Using (4.4.7) we obtain, for the leading order. 



A{q^) oc [< 5p5p > {q^) -q^ < 6p6X > {q'^)] . 



(4.5.3) 




(4.5.4) 




P{R;0,go) 



\AoiR.go)\^ 



NAo{R,go) 



3 1 + 4i?[- cos2 (f ) cos(6'o) + 



cos2(l|I)sin(^o) 
6*0 +7r 



(log(4i?cos2(^)) 



ill)] 



(4.5.6) 



Including the subleading order, we obtain 



P{R-jlgo;Z^)^P{R;0,go)l + 
ARcos^A - cos(^o) - sin(6'o)- 



NAo{R,go)p,lD{Z^,n) 
327r2 

(tt + go)(l + cot(go) tan(^)) - tan(^) "| 2 
l+(7r + ^o)(tan(%)-cot(^o)) J 



(4.5.7) 
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The variable is free and expressible in terms of 6*0 by the second equation in (4.4.7). 
The function D{Z-,^^n) is given in (4.4.12). From (4.5.7) we obtain the cutoff correction 

' 2,^0, , (7r + ^o)(l + cot(^o)tan(|))-tan(|) 

4i<cos (— ) — cos(^o) ~ sm(^o)" 



2 l+(7r + ^o)(tan(f)-cot(6'o)) 

" (4.5.8) 

Its structure is similar to the one seen in the calculation of the width to mass ratio. Note 
that again the dependence on Zf^ and n factorizes and that it comes in through the same 
coefficient as in (4.4.11). Here this is even more non-trivial because it holds for all center of 
mass energies measured by R. This result provides additional evidence for (4.4.13) showing 
that for the cross-section we would be needing the same c as for the width to mass ratio. 
All the dependence on qq and on R is coming from the operator O and the action Sji- 
Another check for eq. (4.4.13) will be provided by our lattice work where we shall see that 
the same dependence on qq and on R enters, the whole difference being expressible by a 
different parameter c. In the present case, a particular observation we can already make 
is that, had we decided to "improve" the width to mass ratio by fine tuning Z^^ so that 
D{Z.j^,n) = 0, we would also have automatically "improved" the scattering cross-section. 

Some feeling for the numbers involved can be obtained from Fig. 7.4d, showing 5|^|2 
as a function of the Higgs mass in TeV for several choices of R and at several values of 

With n < 100 and < Zj^ < 0.5, the cutoff effect on the square of the invariant 
amplitude is less than 3 percent for center of mass energies less than four times the Higgs 
mass, as long as the Higgs mass does not exceed 0.820 TeV. So one can have a strongly 
interacting Higgs particle, at infinite N, with relatively minor cutoff effects. 



4.6. Leading corrections to m/j. 



As promised in subsection (4.1) we present here the leading correction to equation 

(4.1.2) there. We write 

Here to|^q is a function of qjiq of the same form as in (4.1.2). By essentially the same 
methods as before, only this time the coupling we keep constant is gjiQ (as defined in 

(4.1.3) ), we derive a formula for Smuo- It reads 

, 967r2cos(|)C(^7r,^i) 
Smno = D{Z^,n) + ^"^^^ ^' ^ • (4.6.2) 

9R0 
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In addition to completing the calculations in section (4.1) this formula shows us that 
an unphysical quantity will not necessarily have its cutoff dependent corrections propor- 
tional to the function L>(Z7r, n) only. The quantity mjj is unphysical because its definition 
in equation (4.1.1) involved Re[det M^(p^)]. This matrix has a rather complicated 
relation to the correlation functions, and it is obvious that the zero of the real part of its 
determinant is not an "on-shell" quantity. Even if we used fine tuning to "improve" all 
the physical observables by setting D{ZT^,n) to zero the correction to niji would still be 
nonvanishing at leading order in the inverse cutoff. However, without "fine tuning", the 
correction to is quite similar numerically to the corrections to the physical observables. 
Therefore, it made practical sense to discuss mji first, exploiting its relative simplicity, in 
order to get an indication for where in the space of actions the mass bound is likely to be 
larger. 



5. PHASE DIAGRAM FOR THE LATTICE MODELS. 



We are going to investigate the lattice models defined in (2.5.6) for the F4 and hypercu- 
bic lattices and in (2.5.7) for the hypercubic Symanzik improved case. We shall introduce 
a common notation for all three cases with the parameters and symbols having a slightly 
different meaning in each case as explained below. 

Using the constraint = 1 we rewrite the lattice actions, up to additive constants, as 



^ Jx,y °^ Jx 

gx,y is the kinetic term with Fourier transform g(p) = + . . . with the higher order terms 
explicitly given below for each case. To obtain the continuum normalization the fields 
have to be rescaled as = y/r]N{PQ + ^2)^- We therefore restrict ourselves to the region 
Po + p2>0- 

It is convenient to use a lattice type dependent Kronecker delta function, Sx,y = bSx,y, 
where the parameter b is the volume of the Brillouin zone on the particular lattice. On F4 
6=1/2 while on the hypercubic lattice b = 1. 

In the F4 case 77, e, and g are given by 



l2 



Hx)gx,y^y) 



y 



(5.1) 



r/ = 6, e=12, / =2^ 



^ X 



dp 



gx,y — 4:Sx,y - ^ Uy^X 
l=<x,x'> 



B* (27r)2 

1 ^ . (5-2) 
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g{p) was already given explicitly in (2.5.5) and B* is the Brillouin zone of F4 [6]. 

For the hypercubic lattice we are considering two cases: 
The first hypercubic case (HC) is without improvement and has 



77= 2 , e 



= 8, j = Y.^ 



Ip J-T, (27r)^ 



gx,y — ^Sx,y - ^2 ^y,x+n 



(5.3) 



gip) = 2 ^(1 - cos(p^)) =P^-^Y.Pt + O^P^) ■ 
The second hypercubic case (SI) has Symanzik improvement, and the definit 



;ions are 



gx,y — 105a;,y ^ 



'*6 



(5.4) 



v-^ 4 1 
9ip) =2^ ^(1 - cos(p^)) - —(1 - cos(2p^)) 



With the above conventions we can treat all three cases simultaneously. We first 
introduce auxiliary fields in (5.1) to make the dependence on $ bilinear and relax the fixed 
length constraint. 



Si = 7^N{(3o + (32)\ ! ^x)gx,y^y) + r/iV^ / X{x) \ [ ${x)g:,,y^y) 
^ Jx,y ^ J X \.Jy 



+ 



(5.5) 



N 
~2 



ix,y 



— 
2 



with 



,y — V yPo + P2 + 2 i^i^) + ^iy))j 9x,y + p{x)5x,y ■ 
We separate out the zero mode by 



(5.6) 



^(x) = v-\- 



—^vH(x) + —^7t(x), V = - , \v\ = V , V-7f = 0, f H{x) = I 'K^{x) = Q , 
VN VN V Jx Jx 

(5.7) 



and integrate out 7? and H. Then we obtain 

Tr\ogK + v^ J^p-v^ J{p' + lx!g)k-\p' + "IgX!) + ^ (^jX" - 



^2 = T 



(5.8) 
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where A' and p' are, as usual, the non-constant parts of A and p and K now denotes the 
operator from (5.6) restricted to the space of functions whose lattice average vanishes. 



5.1. Saddle point equations and dominating saddles. 



The saddle point equations are 

2 1 /"I ^\^Ps 



viPo + 132 + As) Jp g{p) + Ps viPo + P2 + ^s) 
_2e 1 f gjp) _ 1 . .... ^ ^ ^ 

P2 " iPo + P2 + Xs)Jp9ip)+Ps {P0 + P2 + Xs)^ Ps iKPs)) ■ 

In this equation we have defined 

p = r7(/3o + /32 + A)p . (5.1.2) 

Also, Ji is the lattice integral introduced by [19] for the hypercubic and by [6] for the F4 
lattice. The definition for SI is the same in our notation. By definition, one has J^l = b. 

The symmetric phase is characterized by z;"^ = 0, > and the broken phase by 
v"^ > 0, Ps = 0. A candidate critical point has v"^ = ps = 0. As for the Pauli-Villars 
models the saddle point equations may admit several competing solutions. We shall find 
the dominating saddle by the method used in section (3.1) for the Pauli-Villars models. 
Our presentation will be much more sketchy here. 

Defining 

u = r^{(5Q + (32 + \) , P* = ^, u*=rj{Po + p2) , (5.1.3) 

V P2 



the function 



log[gip) + p] + blogu+{v^ -l)up+^P*{u-u*f . (5.1.4) 




In the symmetric phase, with = 0, the equation d'^ /dp = can be used to define a 
function p{u) for < tt < tto 

Jl{p{u)) = u, wo = <^l(0) = ro , (5.1.5) 

with the lattice constant rg given by 

ro = 0.13823047 for F4 , ro = 0.15493339 for HC and tq = 0.12919024 for SI . 

(5.1.6) 
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The values for the HC and F4 cases can be found in [19,6], while the value for the SI case 
was computed by similar methods. We then derive 

*K,Pa)= / '[G(«)+/3*(«-«*)] + *(«0,0), 

•^"0 (5.1.7) 

Giu) =- - J7^iu) for Q<u<UQ . 
u 

In the broken phase {d^/dp)^—Q = together with > yields the restriction 
Us > uq. Then, extending the range of the function G from u G (0, uq) to the segment 
u e [uq, 00) by G{u) = b/u, we can write the function ^ at the saddle as 

*(«s,0)= / \g{u)+P*{u-u*)]+-^{uo,0) . (5.1.8) 

Again we have identical forms in both phases and can introduce a geometrical inter- 
pretation like in the Pauli-Villars case. We find a tricritical point at 

Pic. = -T = ' ""i.e. = 2«0 = 2ro . (5.1.9) 

«0 ^0 

In terms of the original couplings this corresponds to 



„2 



-r^ , m,t.c. = — 



The second order critical line is described by 



u 



= (f- + l)«0, ^<l, (5.1.11) 



which translates into 



P0,c = + 1) ^ - /32 , P2< P2,t.c. • (5.1.12) 

The restriction {Pq + P2) > leads on the critical line to the requirement 

P2>-P2,t.c. = —^ ■ (5-1-13) 

To ascertain the local stability of our saddle points we need to investigate the small 
fluctuations around the saddle points. This will be done in the next section. 



46 



5.2. Small fluctuations in the broken phase. 



Separating out the zero mode from the field ^ as in (5.7) we want to expand the action 
(5.5) around the saddle point. For this we introduce the shifted fields, recalling that Ps = 
in the broken phase, 

5X = X-Xs, 5H = H, 5tt = tt, Sp = p . (5.2.1) 
We can easily read of the pion propagator in Fourier space 

< (57r"(57r^ >= — — ^ . (5.2.2) 

r}{l3Q + 132 + Xs)g{p) ^ ' 

From this we find the pion wave function renormalization constant 

We see that the renormalized pion propagator is independent of the couplings and hence 
so are its cutoff corrections. We also see that in the SI case the inverse pion propagator 
has no order contribution and thus is Symanzik improved. 

Using given above and — Zj^f^ we obtain from the first saddle point equation 
in (5.1.1) {ps = in the broken phase) 

^^ = ^72 • (^•2-4) 
ro + f4 

Then, also using the second saddle point equation, we find 

Therefore, at fixed we can trade the coupling Pq for 

Integrating out the pions from the action (5.5), with the zero mode separated out, we 
obtain, up to quadratic order in the fluctuations. 



Si^^ = \ I ^aMab^^B ■ (5.2.6) 



P 



Here the three component fields V'A represent 5H^ 5p and 5\ for A= 1,2,3 respectively. 
The non- vanishing entries of the symmetric matrix Mj^b(p) are, neglecting 1/N correc- 
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tions, 

Mi2 = "//Vv 
^13 = 



1 " -v2 



M22 = - -iV^^ / ^ = — NZillp) 

M23 = - -NZtr] / ^^^^^ ^ — ^ , ^ = — NZ^rjro 

2 Jk g{k+ ^p)g{k- ^p) P2 2 V ^^^7^P2/ 

(5.2.7) 

We need two "bubble" integrals, I{p) and (5(p)- -^(p) has already been computed in an 
expansion in the momentum p in [19] for the HC lattice and in [6] for F4. The computation 
for SI, employing the same method, is straightforward. One finds for F4 and SI 

1 



I{p) = -j^ logp^ + ci - C2p2 + 0(/logp2) (5.2.8) 



(5.2.9) 



with the constants q 

ci = 0.0466316 , C2 = -5.4968 ■ 10"^ for F4 
ci = 0.0283716 , C2 = -3.5981 • 10"^ for SI . 

For HC the result is 

Hp) = -j^log^'^ + - C2,ip^ - C2,2 J^P^/P^ + 0(/log:p2) (5 2.10) 

with 

ci = 0.0366783 , C2,i = -7.524 • 10~^ , C2,2 = -2.6386 • 10~^ . (5.2.11) 

Again, at order 1/A^, the HC lattice has a Lorentz invariance breaking term. However 
we shall need I{p) only for 'time' like, and thus on-axis, momenta. Then the HC result 
reduces to the form (5.2.8) with C2 — C2,i + C2,2 = —3.3910 • 10~^. We see here explicitly 
how the specific choices of observables that we made effectively hide the Lorentz breaking 
effects of the HC case. 

An expansion in p for the integral Q{p) is easily computed with the result 



Q(p) = 6 + 7p2 + 0(/) 
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(5.2.12) 



with 



7 = 0.00661524 for F4 , 7 = 0.01496670 for HC and 7 = 0.01736345 for SI . 

(5.2.13) 

We are now in the position to discuss the local stability of our saddle points. It is easy 
to see that on the critical line, i.e., for Z^^ = Tq-^ ((5.2.4)), M33 vanishes ai p — when (32 
approaches P2,t.c. from below. Atp — Mn and M13 also vanish and hence at the tricritical 
point det M = 0. Thus, as in the Pauli-Villars models, the Higgs particle becomes massless 
at the tricritical point and also plays the role of a dilaton. For (32 < P2,t.c. we find, at least 
in the neighborhood of the critical line, that det M stays positive for Euclidean momenta. 
Thus the saddle points in the critical regime that we are interested in are locally stable. 



6. HIGGS MASS BOUND AND CUTOFF 
EFFECTS WITH LATTICE REGULARIZATIONS. 



As in the Pauli-Villars case we will first do an approximate calculation to see the basic 
trends. Next we follow up with an accurate evaluation including the computation of the 
leading cutoff correction to the width to mass ratio and to tt — tt scattering. 



6.1. Approximate calculation. 



We will first study the quantity m^, defined by 

Re[detM(p)]^^(^^^,^^ = . (6.1.1) 

From (5.2.7) we obtain 

detM(p) = jN^Z^r]'^ 

2. \ / X 1 1 (6.1.2) 

Qip) - ^2^J [9{p)m + Vl) - 9{p)rl - 2rof^g{p) + -flg{p?m • 

We use Z~^ = rQ + and solve (6.1.1) for m|j,, << 1, counting logs to be of order 
1. To leading order we find, as expected, the form (4.1.2,4.1.3) with the proportionality 
constant C now depending on P2 

C(/?2) = exp ^ Stt^ci ^} . (6.1.3) 



VW2 
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We see that C diverges, as for the Pauh-ViUars models, when the tricritical point is 
approached. Again we want to be as far away from this point as possible. 

When we decrease /32, 0(^2) and hence mji (in lattice units) becomes smaller at 
constant gji (which means that the mass is kept fixed in physical units although mji 
varies). Then the cutoff effects become smaller and gji can be allowed to increase. Thus 
we expect to obtain the largest mass bound at the smallest acceptable P21 (5.1.13). The 
ratio of the C's for the standard non-linear action and for this extreme case is found to be 

^ ^ = exp <^ — — ^ } . (6.1.4) 



C{-P2,t.c.) 

This evaluates to 4.521, 2.580 and 1.933 for F4, HC and SI respectively. According to 
the rough estimates of section (4.2) we would expect the change of /32 from to — /92,tc to 
induce increases of about 0.075 TeV, 0.047 TeV, 0.033 TeV in the Higgs mass bound for 
the three cases. In the next section we shall see that these rough estimates for the increase 
of the Higgs mass are accurate within a factor of two. 

As in section (4.6) for the Pauli-Villars case we can compute the leading correction to 
mji given by (4.1.2) and (6.1.3). In the notation of (4.6.2) we find 



where C, is defined through g{p) = — (^p^ + 0{p^) for on-axis momenta p. From equations 
(5.2, 5.3, 5.4) we find C = 1/12 for F4 and HC and C = for SI. This correction is less than 
15% in the region of interest. 

We can make the lattice result look even more like the Pauli-Villars one by going to 
continuum normalization. We denote the corresponding pion wavefunction renormalization 
constant by Z.j^, Z.j^ — r]{(3Q + ^2)^^- Then, using eq. (5.2.4) and (5.2.5) we obtain 

V{P0 + P2) = Z^{ro + fi) 

b P2 ^ ^TT-l (6-1-6) 
2e (/?o + P2? ~Zl 

Therefore we can now trade the couplings /3o and /32 for and Z^r- The critical line is 

now given by = and < < 2. The upper limit corresponds to the tricritical 
point. = 1 denotes the standard non-linear lattice action, without the term with four 
derivatives and four fields. 

In this normalization the term Q{j)) — ^2^2^^ in det M of cq. (6.1.2) has to be replaced 

(Z —2) 

+ Q{p) ~ b. The proportionality constant in eq (4.1.2) now depends on and 
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becomes 

C{Z^) = exp { Stt^ci - } (6.1.7) 



b - 2 

to be compared to eq. (4.1.4) for the Pauli-Villars models. 

However, keeping fixed, while varying does not correspond to keeping /?2 fixed. 
In a numerical simulation of the lattice models one cannot easily follow a line of fixed 
Ztt. It is much more natural to carry out simulations at a fixed value of /?2, approach the 
critical point along this line and subsequently vary /?2. Therefore we will return in what 
follows to the usual lattice normalization. 



6.2. Leading correction to the width to mass ratio. 



Again, the computation in this section closely follows that of the Pauli-Villars model, 
section (4.4). The Higgs resonance is a complex root on the second sheet of detM(p) = 
for p = {i{mjj — ^7if)) 0)- det M is given in (6.1.2) where we again insert = tq + f^. 
The leading order solution, with notation (4.4.5) and (4.4.6) is 



2tan(6'o), 
967r2cos2(^)sin(^o) 



1^0 = CW exp 



90 



(6.2.1) 



with C{(32) given in eq. (6.1.3). 

To compute the sub leading terms we expand again as in (4.4.9) and set 6g — (see 
(4.4.10) and the explanation leading to it). Straightforward algebra then leads to 



167r^ 



S {go; P2) 



(^2 (6,;2^2-2er2^2 



^0 



M 



H 



1 + (tt + ^0) [tan (^^) - cot(^o) 



(6.2.2) 



As for the Pauli-Villars case, the cutoff correction factorizes into a universal g depen- 
dent part and a part carrying all the information about the cutoff through the dependence 



on /?2. The latter is given here by 167r [c2 



C (/32). The universal factor is 



identical to the one obtained in the Pauli-Villars case, eq. (4.4.11). 

On the lattice we are also interested in comparing directly to numerical data obtained 
at = 4. For this comparison we need some "unphysical" quantity, like mjj. We therefore 
write down the explicit formula to first subleading order in the cutoff 
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(6.2.3) 



and obtain 



™iJO = /^O cos 



IGtt C2 



i4v Pi 



COS ^0 + sin(^o) 



(tt + 6*0) (1 + cot(^o) tan(%)) - tan(%) 



(6.2.4) 



l + (7r + 6'o)(tan(f ) 



cot(^o)) 



+ 



sin ^0 



c- 



broV Po 



2er2)2 



6.3. Leading correction to tt — tt scattering. 



Scattering on the lattice was discussed in [19] for the HC lattice and generalized in [6] 
to the F4 lattice. It is easy to see that at first non- vanishing order in 1/N the invariant 
amplitude A{s,i,u), where s, i and u are the lattice versions of the Mandelstam variables 
s, t and u, [6] describing tt — tt scattering in the center of mass frame, is only a function 
of the lattice center of mass energy squared s. For our purposes we do not need the full 
expressions for the variables with hats and it suffices to know that s = W"^ + 0{W^) where 
W is the center of mass energy. The differential cross-section with equal isospin indices on 
the incoming pions is given by 



da\ N 



■\A{s)\^ . (6.3.1) 



The kinematic prefactor has no cutoff corrections of order l/A^. For the F4 and SI models 
this is true for incoming and outgoing momenta in any direction, while for the HC model, 
due to the Lorentz invariance breaking at order 1/A^, this property only holds for on-axis 
momenta. For the F4 and SI models anisotropics in the differential cross-section only occur 
at order 1/A^. 

To evaluate A{s) we need to know the pion interactions. For this we go back to the 
action (5.5) and introduce pion and Higgs fields via eq. (5.7). We then expand around the 
saddle but do not integrate over 5H and 5tt yet. The part of the action giving the pion 
interaction is then 



dTv{x)gx^yS7v{y) 



+ i / 5pix){5nf{x) . (6.3.2) 



The TTTTp and ttttA vertices are now easily obtained, and one can see that the latter vanishes 
for on-shell pions, since then g{p) = 0. Thus, the on-shell amplitude with q = Pi + P2 is 
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given by 



A{q) oc < 6p6p > (g) 



(6.3.3) 



Following closely section (4.5) where the Pauli-Villars case was dealt with and using the 
same notation, we find the same leading order result for the invariant scattering amplitude, 
eq. (4.5.6). Including the subleading order we obtain on the lattice 



= C2 



(6772/32 - 2er2)2 



l4Re[NAo{R,go)] 



4i?cos2(|) 



cos(^o) ~ sin(^o) 



(7r + ^o)(l + cot(6'o)tan(f )) 



tan(f) 



l+(7r + ^o)(tan(f)-cot(^o)) 

(6.3.4) 

This cutoff correction has again factorized into a cutoff dependent part and a function of 
qq which is the same as the one found for the Pauli-Villars models, eq. (4.5.7). The cutoff 
dependent part is identical to the cutoff dependent factor we have obtained for the width 
to mass ratio in equation (6.2.2). 



7. FROM AT = oo TO AT = 4. 



In this section we shall extract quantitative information about the Higgs mass bound 
from our large results. To do that we need to first discuss how large the corrections to 
the N = oo limit are when N = A. Our objective will be to argue that, in the region of 
interest for the bound, the N = oo numbers should be expected to be different from the 
N = A values by about 25 percent (l/N). The question of the difference between N = oo 
and N = A has been studied in some detail before by Lin, Kuti and Shen[20]. They arrived 
at the conclusion that the difference between N = oo and N = A is very large and that 
N = oo results are qualitatively wrong at N = A. We disagree with this conclusion in the 
region where the bound is close to saturation, and this is the region of interest to us. Some 
of the approximations used in [20] for the lattice "bubble" diagram are unjustified and we 
suspect this to be the reason for the discrepancy between our views on the usefulness of 
the expansion for the Higgs mass bound problem. 

Throughout this and the next section we undo the rescalings of F-,^ and f-jr by VN so 
that, at iV = 4 = 0.246 TeV. 



7.1. Simple error estimates. 



The quantity that we are extracting quantitative information from is the tt-tt scattering 
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invariant amplitude ^(9 )• From it we extracted the Higgs mass and width, mjj and 7^, 
the quantity mji and the cutoff correction ^|^|2. Moreover, from the behavior of the 
amplitude when q'^ ^ wc can also extract f^. Therefore our basic question is how 
accurately is the "true", = 4, scattering amplitude A{q^) approximated by the N = 00 
expression for complex with < 4m|^. 

The best approach would be to compute the 1/N correction explicitly but this is a 
demanding calculation which we have not done. We are attempting to guess what the 
right order of magnitude of the result of such a calculation would be. The most naive 
guess is that the leading correction will be suppressed by 1/N, meaning that our results 
are good to 25% when applied to = 4. Indeed, the most conspicuous omission in the 
N = 00 limit is the fact that the number of pion types is — 1 rather than A^; the 
naive error estimate then follows from the observation that pairs of pions of identical type 
contribute more or less additively. There is one example where this sort of pion counting 
argument is believed to be exact: When one computes the finite volume leading correction 
to the < (f) ■ (f) > correlation at zero momentum in the 1/N expansion, the leading and 
subleading terms are related by exactly a factor of N [21]. Moreover, it is believed that all 
higher order corrections in 1/N vanish if the pion wave function renormalization constant 
is extracted [22] . 

It is well known [23] that one cannot always expect the 1/A^ corrections to obey these 
naive estimates; this point was also discussed in ref. [20]. Let us repeat the argument here, 
but rephrased for the amplitude A{q'^): If the system is sufficiently close to criticality, the 
low q"^ behavior of A can be extracted from an appropriate RG equation and is dominated 
by the noninteracting fixed point. The function P{g) appearing in this equation can be 
approximated by its leading term in g and that term is calculable because it comes from one 
loop diagrams. It is known that, to one loop, the beta function for the properly rescaled 
coupling is proportional to 1 + 8/N and, thus, the finite N answer is for A?" = 4 larger by 
a factor of 3 than the N — 00 answer. This will affect the amplitude because the solution 
of the RG equation will simply say that A{q'^) is well approximated by the tree level 
expression with an effective coupling that is dependent on q^ with a dependence dictated 
by the /3-function. If the coupling is also weak at cutoff scales, the effective coupling can be 
related to the bare coupling via the one loop P function; in this way one obtains A{q'^) in 
terms of the bare coupling. If we are very close to the critical point the cutoff is very large 
(for example, when compared to the Higgs mass) and the effective coupling becomes almost 
entirely independent of the bare coupling; in that case the large error in the /? function 
induces a large error in A(q'^) and we can't trust the numbers obtained ai N = 00 when 
N = A. Note that the difficulty arises due to the existence of two very disparate scales in 
the problem: One scale is the cutoff and the other is set by the pion decay constant. Two 
very different scales exist only when the system is very close to a critical point because /tt 
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vanishes there. 

For the Higgs mass bound problem we are mostly interested in cases where ^ttFtj- A; 
in these cases the usefulness of the RG equations is very limited because we do not have to 
connect two very different scales. Therefore, for those actions that lead to the approximate 
saturation of the bound we do not need to worry about the 8/N error in the /3-function 
and can hope for accuracies of order 25% at N = 4. In the next section we shall use 
available numerical data as evidence that this hope is realized in the broken phase of the 
models we are interested in. 

In the symmetric phase more is known about the 1/N series. For the simplest nearest 
neighbor action on the hypercubic lattice one is able to get into the region where the 
triviality bound on the self-coupling is saturated (at a correlation length ^ of about 2 
lattice spacings) with an accuracy on the coupling of the order of 2-3% if one keeps three 
orders in at = 4 [24]. This is compatible with the assumption that the 1/N series 
for the coupling has a radius of convergence of order unity when ^ ~ 2. 

Also, regarding the value of the bound on the coupling in the symmetric and broken 
phases, one can compare the results obtained for the single component model (A^ = 1) to 
those obtained for = 4; in spite of the fact that there are no Goldstone bosons at = 1 
the upper bound on the coupling is relatively stable when expressed as a fraction of the 
tree level unitarity bound [25,26,19]. Again the 1/N series looks reasonably well behaved. 

7.2. Numerical test of error estimates. 

In Figures 7.1a to 7.1c we compare the large N predictions with A^ = 4 results for 
the simplest lattice models investigated to date by Monte Carlo and independent non- 
perturbative means. For the F4 lattice (Fig. 7.1a) the numerical data follows the large A^ 
prediction quite consistently over the entire range of interest, 0.3 < rrifj < 0.8, with the 
large A^ numbers exceeding the A^ = 4 numbers by 20-30%. For the hypercubic lattice the 
excess is slightly smaller, 15-20% in the region of interest, and for the hypercubic improved 
less than 15%. 

More importantly we see that the A^ = 00 difference in Mjj/F-jy between the various 
lattices in the region of interest underestimates the A^ = 4 difference by up to a factor of 
two. Therefore, if we knew the A^ = 4 results for only one of the lattices we could predict 
the results for the others, using infinite A?", to an overall relative accuracy of about 7-9%, 
which is of the same order of magnitude as the statistical error in a typical Monte-Carlo 
simulation. One of our objectives is to predict the results for the coupling at AT = 4 for new 
lattice actions. Being interested in actions that are not fundamentally different from the 
ones that have already been investigated by other means, we can expect our predictions 
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Figure 7.1 rnjj/fT^ vs. mjj for the naive {^2 = 0) lattice actions. The sohd hne is 
the large N result scaled to N = A. The diamonds in fig. (a) are numerical 
results from [27]. In fig. (b) the boxes are results from [19], the diamonds are 
numerical results from [28] , and the crosses from [29] . In fig (c) the diamonds 
are numerical results from [29]. 

based on large N to have an accuracy of the order of several percent. Preliminary results 
in the F4 case provide further support to this claim [2]. This is very useful because one 
does not want to invest resources in simulating actions that end up not giving a higher 
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bound. We believe that quite reliable order of magnitude estimates (and definitely the 
right sign) for the effects of different terms in the action on the bound can be obtained 
from the N = oo results. 

Until now we only tested the large N predictions against results that could be obtained 
at A?" = 4 by other non-perturbative means. But the estimates for the bound also need non- 
perturbative evaluations of the cutoff effects contributing to physical observables. There do 
not exist any tested, practical, non-perturbative methods yet that make it feasible to obtain 
numerical evaluations of physically observable cutoff effects at = 4. ft is possible that 
such methods will be developed but it is doubtful that one will be able to reach reasonable 
accuracies for estimating such small corrections by Monte Carlo methods. Therefore we 
are left with only two options for evaluating the cutoff effects: the loop expansion and 
the 1/N expansion. It is gratifying that the order of magnitude of the results at tree 
level order for nearest neighbor hypercubic and F4 actions are in agreement with the large 
N estimates when the coupling constants are close to the upper bound. We could have 
worried about the tree level estimates because, as explained in Appendix B, higher loop 
effects are not negligible and a method for resumming the "leading log" terms contributing 
to order terms is not yet available. However, it seems that this does not have a major 
effect when the coupling is sufficiently large, probably for the same physical reasons that 
the /9-function problem was relatively harmless for large couplings. We are further helped 
by the fact that the bound depends weakly on the magnitude of the cutoff effects: Changes 
of one order of magnitude in the latter have an effect on the bound of the order of only 
a few percent. If we assume that, when the coupling is sufficiently strong, the accuracy 
of the N — 00 answer is of order 25% for the cutoff effects too, we can use the large N 
estimates for the leading cutoff effects to estimate the Higgs mass bound in a Monte Carlo 
calculation. We have not devised a method for computing the cutoff effects in the loop 
expansion for the more complicated actions that were the subject of this paper because we 
feel that the large N results are at least as reliable. Nevertheless, we think that it would be 
interesting to see the explicit N dependence emerging from a perturbative computation. 



7.3. The large N numbers. 



Let us now go over the quantitative results of our work. We start with the kind of 
results that are the typical outcome of a Monte Carlo simulation. In Figures 7.2a-d we 
show plots of the ratio mj^/ f-^- as a function of the quantity mjj. The latter is equal to 
Mjj/A and contains an implicit understanding for what the cutoff is taken to be. For all the 
lattice results we take the cutoff as the inverse lattice spacing of the elementary hypercubic 
lattice on which the action is defined (the F4 lattice is viewed here as a hypercubic lattice 
for which a choice of action was made so that no fields reside on sites who have an odd sum 
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of integer coordinates). In the Pauli-Villars case the cutoff was defined as the absolute 
magnitude of the distance, in the complex momentum square plane, to the ghost pole in 
the pion propagator closest to the origin. 
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Figure 7.2 mjj/ f-j^ vs. mij: The solid lines represent the naive actions (/?2 = on the 
lattice and Z-,^ = 1 for PV) the dotted lines the actions with a four derivative 
term turned on to maximal allowed strength (/?2 = —P2,t.c. on the lattice 
and = for PV). In figures (a)-(c) rriff is measured in lattice units and 
in fig. (d) it is measured in units of Ag. 
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For the three lattice cases we show two hnes on each graph. The lower one gives 
mji/ Jti for the simplest action which is bilinear in the fields. From the higher line one 
can read off the increase in mjj/ fn that can be induced, at the same value of Mjj/A, by 
turning on the four field terms to maximal acceptable strength in the direction of increasing 
repulsion between the pions. We see that the largest relative effect takes place in the F4 
case while the smallest relative effect takes place in the hypercubic "Symanzik tree level 
improved" case. Since the lines for the simple actions give higher values of mjj/ for 
exactly those cases for which the effect of the four field terms is smaller, we see a tendency 
towards "equalization" when the four-field terms are turned on. This fact supports some 
sort of approximate "universality" among the various bounds. Of course we cannot expect 
anything rigorous to hold in this respect, but it is impressive that the higher lines for all 
three lattice actions are closer to each other than the lower lines are. 

In the Pauli-Villars case we can compare different values of n and see similar effects. 
While the overall "plateau" attained by the higher curves is quite similar to the lattice 
detailed comparison cannot be made because there is no "natural" relationship 
between the "cutoffs" in the two schemes.* 

There is a simple explanation for the systematics we observed: Suppose each one of the 
actions is evaluated for slowly varying fields and a field rescaling is carried out in accordance 
with eq. (2.5.2) to remove the four derivative term from the term in the Lagrangian that 
is bilinear in the fields (this cannot be done for the simple hypercubic action and we shall 
discuss this presently). We see that the ordering of the lattices by rnu/ f-,^ at fixed mjj 
and with naiVe actions follows the magnitude of the four field coupling induced by the 
field redefinition of (2.5.2). For the hypercubic lattice there is no field redefinition that 
can eliminate the four derivative term but it is clear that this term is "weaker" than the 
corresponding term on the F4 lattice, and obviously "stronger" than the absent term in 
the improved case. We obtain for the first time an explanation for why F4 results for 
the bound have turned out to be smaller than the hypercubic ones [27] and preliminary 
numerical results with "Symanzik improvement" seem to give bounds [30] that are larger 
than the ones obtained on hypercubic lattices with the simplest action.^ 

Similar logic works when one tries to understand the effects of the four field terms. 

* At = cxo one can define the ratio of the Pauli-Villars cutoff to a particular lattice cutoff 
by requiring exact matching of the coefficient c in equation (4.4.13). This is equivalent to 

an exact match of the cutoff effects at order 1/A^. However, this match cannot always 
be realized since some Pauli-Villars actions can make c vanish and this is impossible to 
achieve with the lattice actions that we investigated. In the text wc meant by "natural" 
the sort of relationship one would guess by just looking at the bare action, not the more 
sophisticated relation described above in this footnote. 

^ Very recent work [29] reaffirms the trend; the newer data points have been included in 
Figures 7.1. 
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The bottom line is that one has full support for the "rule" that strengthening the repulsion 
between the pions in the low momentum approximation to the action, viewed as a chiral 
effective Lagrangian, increases the bound. Still, the answer is not entirely universal, and 
there are some finer differences, especially between the Pauli-Villars cases and the lattice 
cases. 

We now turn to the cutoff effects. In order for the features observed above to really 
imply that higher Higgs masses are possible when the action is changed, we have to show 
that the above trends are also obeyed by the cutoff effects. This would free us from the 
somewhat arbitrary choice of comparing different lattice actions at the same mjj with 
a specific choice for the "cutoff" in each case. Because the mechanisms explaining the 
observed trends, even with this arbitrariness unresolved, are of a physical origin it is to 
be expected that the cutoff corrections will indeed follow these trends. Therefore just for 
the purpose of knowing which kind of action one should try to simulate at N = A the 
indications obtained up to this point should be taken seriously. However, to make well 
defined statements we shall need to know the magnitude of the cutoff effects in greater 
detail anyhow, even after Monte Carlo data replaces some of the large N graphs in Figures 
7.2a-c. We therefore show in Figures 7.3a-d the cutoff effects on the width to mass ratio 
for the four cases under consideration. The numerical value of these effects is rather small, 
and this is in agreement with tree level perturbation theory in the simplest lattice action 
cases. We prefer to use the more stringent cutoff effects on the differential cross section for 
TT-TT scattering at ninety degrees in the center of mass frame. We chose this observable to 
be able to compare with existing perturbative results. The cutoff effects on the invariant 
amplitude square are shown in Figures 7.4a-d. Note that the trends are the same as in Fig 
7.3a-d. 

We see that the trends already seen in Figures 7.2a-d are respected. The amount of 
possible increase in the bound is smallest for the hypercubic improved case and largest for 
the F4 case. We can immediately read off from the horizontal axis the mass in TeV that 
is permissible for the Higgs mass if one puts some restriction on the cutoff effects. Again 
one sees a tendency toward approximate "universality" with a bound on the Higgs mass at 
N — 00 oi about 0.820 TeV* This is achieved with the smallest cutoff effects in a Pauli- 
Villars type of regularization. What seems to be special about our PV regularizations is 
that the pion propagator has no order p'^ term; this is also true of the SI action and can 
also be achieved on the F4 lattice, but at the cost of introducing next nearest neighbor 
interactions. t The following is an "improved" action on an F4 lattice with a tunable 

* The figure also shows that one does not need to set P2 exactly —(^2,t.c. to obtain an 
effect of similar magnitude. 

^ While it is true that order terms in the bare action can be removed by field redefinitions, 
the nonlinear relationship between the bare parameters in the action and the cutoff effects 
implies that the elimination of the terms by field redefinition is not exact. 
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Figure 7.3 Leading order cutoff effects in the width to mass ratio: The sohd hne repre- 
sents the naive actions (/?2 = on the lattice and Z-,^ = \ for PV) and the 
dotted hne the actions with a four derivative term turned on to maximal al- 
lowed strength (/32 = — /?2,i.c. on the lattice and Z-^- = for PV). The dashed 
line in the PV case shows an example with the "wrong" sign in front of the 
four derivative term, corresponding to Z-j^ = 2.5. 
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Figure 7.4 Leading order cutoff effects in the invariant tt — tt scattering amplitude at 90^ 
for the naive actions (/?2 = on the lattice and Z^^ = \ for PV) and for the 
actions with a four derivative term turned on to maximal allowed strength 
(/32 = —P2,t.c. on the lattice and = for PV). The dotted line represents 
center of mass energies W = 2Mjj, the dashed line W = SM^^ and the solid 
line W = AM In the PV case we again show an example with the "wrong" 
sign in front of the four derivative term corresponding to Z-,^ = 2.5. 
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"triangle term" 



S= -2N{(3o + pi + P2) 



2 J2 Hx)-i{x')-^ Yl Hx)-Hx'") 
<x,x'> 

-NPi Yl [^(^) • ^(^') - 1]^ 



<x,x"'> 



<x,x'> 



8 



E 



x <U'> 

ina;V0, i'na;'V0> x,x',x" all n.n. 



^{x) ■ ^{x ) - l] U{x) ■ ^{x") - 1 



(7.3.1) 

where "' denote sites, < x, x' >, Z, l' links, and < x, x'" > next nearest neighbor- 

ing pairs. The field is constrained by ^^{x) = 1. In (7.3.1) only sites separated at most a 
distance of two, in units of the embedding hypercubic lattice, are coupled. 




-0.03 -0.02 -0.01 0.00 

Figure 7.5 Leading order cutoff effects in the invariant 90^ tt — tt scattering amplitude 
vs. /?2 for three values of Mjj on the F4 lattice. The center of mass energy is 
set atW = 2Mh. The cutoff effects for a 0.815 TeV Higgs at P2 = -p2,t.c. 
and for a 0.720 TeV Higgs at /?2 = are equal in magnitude. 
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To see more clearly how the bound changes with the action we show an example in 
Figure 7.5 for the F4 lattice: identical cutoff effects, of 4%, will be found at AT = 00 for 
a Higgs mass of 0.720 TeV with the simplest action, and for a Higgs mass of 0.815 TeV 
with an action that has the maximal amount of four field interaction we have allowed. It is 
important to understand that this difference is substantial when the width is considered. 
In Figure 7.6 we show the large N Higgs width as a function of the Higgs mass (note that 
the leading order weak coupling approximation severely underestimates the width when 
the Higgs becomes heavier). In our example the width went from 0.320 TeV to about 
0.500 TeV and the heavier Higgs is definitely strongly interacting. There is no doubt 
therefore that, at least at infinite A^, stopping the search for the Higgs mass bound at the 
study of the simplest possible actions would have been misleading. 
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Figure 7.6 The regularization independent part of the width vs. Mjj. The solid line 
displays the large N result scaled to N = A and the dotted line shows the 
leading order term in perturbation theory. For large Mjj, the perturbative 
answer severely underestimates the width. 
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7.4. Predictions for N = 4. 



We now turn to the most speculative, but quite important part of our work, namely 
numerical predictions for the physically relevant case at N = A. Our approach is as follows: 
For the simplest actions (which are bilinear in the fields) numerical results at N = A are 
already available. We compare the Mjj/Fjj- ratios we calculated at A?" = oo, viewed as 
functions of mjj, to these numbers. We then use the difference between the N = oo and 
N = 4 numbers, at a given mjj, as an estimate for the same difference for each one of the 
new actions that have not yet been studied numerically. In this way we are able to make 
predictions for the new actions at = 4 based on our N = oo results. We can test this 
method by assuming that we have numerical results for only one of the simple actions and 
"predict" the numerical results for another simple action. We find that our "predictions" 
are good to a few percent in the region where cutoff effects on the tt — tt scattering cross 
section are of the order of a few percent. Comparisons between the N = oo and N = A 
results for the simplest lattice actions are given in figures 7.1a for F4, 7.1b for HC, and 
7.1c for SI.* 

For the F4 lattice, we see from figure 7.1a that for the range of interest, ^ mjj > 0.5 
(with rrifj the Higgs mass in lattice units of the large N calculation at N = A), the large N 
results are larger than the numerical ones by less than 0.150 TeV. For the HC lattice, figure 
7.1b indicates that for mjj > 0.4 the large N results are larger by less than 0.100 TeV. 
Finally, for the SI case, figure 7.1c shows that for rriff > 0.3 the large N results are larger 
by less than 0.110 TeV. From these figures we also see that the large results for are 
consistently larger than the numerical ones. Therefore, given a large N estimate for Mjj, 
we predict that the corresponding value, calculated in a numerical simulation at = 4, 
will not exceed the large A^ estimate. We also predict that it will not be lower by more 
than the amounts given at the beginning of this paragraph, provided that mjj stays within 
the ranges of interest indicated there. 

The largest Higgs mass is obtained when the four derivative term is turned on to 
maximal allowed strength. Using the guidelines of the previous paragraph and figures 7.3a 
- 7.3c or 7.4a - 7.4c we can estimate the Higgs mass bound for a given amount of allowed 
cutoff effects. We decide to allow up to 4% cutoff effects in the cross section at W < 2Mjj 
as computed in large A^. For the F4 case we find that the bound should be between 0.660 
and 0.810 TeV; this prediction is in agreement with our preliminary data [2]. For the HC 
it should be between 0.740 and 0.840 TeV, and for the SI between 0.750 and 0.860 TeV. 

* Our plots include very recent data that appeared in [29]; this preprint became available 
during the final stages of preparing our paper for submittal. We have not yet had the time 
for a thorough study of [29] . 

^ This is the range in which the bounds are obtained when —P2,t.c. ^ P2 
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The second number in each case is the large N result ignoring the likely overestimates and 
hence is a very conservative estimate for the bound. 

Although we do not have any large N results corresponding to the action (7.3.1), we 
expect the bound to increase by about 0.050 TeV from the value given for the F4 lattice 
in the previous paragraph, to somewhere between 0.710 and 0.860 TeV. This is expected 
because the action (7.3.1) should be closer to the PV case, and from figures 7.4a and 7.4d 
we see that the Higgs mass bound given by the PV case is larger than the one given by 
the -F4 lattice by at least 0.050 TeV. 

We see that if we allow up to 4% scaling violations in the cross section atW< 2Mff, we 
should expect, without being overly conservative, a Higgs mass bound of about 0.750 TeV. 
Renormalized perturbation theory suggests that such a heavy Higgs particle will have a 
width of at least 0.210 TeV. The large N results indicate an even larger width, around 
0.290 TeV, already subtracting a possible overestimation by about 25%. It is unclear 
whether such a wide Higgs particle can still be seen on the lattice without applying more 
sophisticated techniques than mere direct measurement. 

The approach of [29] is not based on direct measurement of the mass in the / = 0, J = 
channel but relies quite heavily on the applicability of ordinary perturbation theory. We 
know from our large N computations that ordinary perturbation theory does not work 
that well for Tjj/Mjj when Mjj/Ft^ gets close to the bound. However it is not ruled 
out that the particular quantities that are calculated in perturbation theory in [29] do 
have a "better" perturbative expansion. Using the expression for the variance of the 
magnetization squared in [4]"^ (evaluated there also for the purpose of finding alternative 
definitions for the scalar self-coupling, definitions that lead to quantities that are more 
easily accessible numerically and don't suffer from finite width contamination) one could 
test the validity of perturbation theory in the precise setting of [29] , at least at infinite N . 

It is important to stress that all our predictions in this subsection are relative to the 
presently generally accepted numerical values at = 4. Thus, if there are systematic 
errors in the iV = 4 results, for example due to ignoring finite width effects in the direct 
Monte Carlo measurements and/or due to the application of perturbation theory in [29], 
these errors will propagate into our predictions. 



8. SUMMARY. 

In this paper we have studied the regularization scheme dependence of the Higgs mass 
triviality bound in an 0{N) scalar field theory to leading order in 1/N. All possible 

^ Equations (6) and (7a-b) there; see also the related discussion. 
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leading cutoff effects can be induced by adding only a few higher dimensional operators 
with adjustable coefficients to the standard action. We have found that the highest 
Higgs masses are obtained in the non linear limit, at infinite A. This has led us to consider 
non linear actions with the leading cutoff effects induced by four derivative terms with 
tunable couplings. 

We used one class of continuum regularization schemes, of Pauli Villars type (PV), 
and three kinds of lattice regularizations: F4, Hypercubic (HC), and Symanzik Improved 
Hypercubic (SI). For these regularization schemes the action, when expanded for slowly 
varying fields to order momentum to the fourth power, is of the form 

(2.5.1) 

where ^ = Nj3. There are four control parameters in this action. One is redundant since 
to this order it can be absorbed into the other parameters by a field redefinition: The 
parameter 60 can be absorbed in hi and 62 by 

$ - ^^=tM!L== . (2.5.2) 



However, eliminating the dependence in 60 to order momentum to the fourth power in the 
bare action does not necessarily imply that the dependence in 60 has been eliminated to 
leading order in the inverse cutoff from physical observables. The vacuum fluctuations are 
"aware" of the full bare action and will carry that information to the non universal part of 
the physical observables. For example the parameter n of the PV case, although associated 
with an operator of dimension larger than four, does appear in the relation between bare 
and renormalizcd parameters and in the leading cutoff corrections to the width and cross 
section. Similarly for the lattice regularizations, some dependence on the full structure of 
the lattice propagators comes in through the constants rg, ci, C2, and 7. Because of this, 
the effects of 69, after the field redefinition, are not completely absent from the leading 
cutoff corrections to the physical observables. Still, the effect of 60 is probably small and, 
as far as the value of the Higgs mass bound is concerned, one may be able to cover the 
whole range of leading cutoff effects by varying hi only. Our results indicate that this 
is realized to a good extend, leading to an approximate universality of the Higgs mass 
triviality bound. 

For PV we simply set 60 = 0) ^^nd for the lattice regularizations 60 is set to the 
naive value obtained in the expansion of the lattice kinetic energy term. We calculated 
the phase diagram for each regularization and found a second order line that ends at a 
tricritical point where a first order line begins. We studied the physically interesting region 
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close to the second order line in the broken phase. The parameter /? corresponds to the 
relevant direction and is traded, as usual, for the pion decay constant /t^. The parameters 
bi and 62 control the size of the leading order cutoff effects for a given value of /tt- Our 
analysis has shown that, to this order, the quantities we considered do not depend on the 
parameter 62- Therefore, a very simple situation emerged with the scale set by /jr and the 
leading cutoff effects parametrized by only one parameter bi. This simplification at infinite 
N makes it easy to relate very different regularization schemes, and leads to a reasonably 
"universal" bound on the renormalized charge g. 

We calculated the leading cutoff effects on two physical quantities, namely the width 
of the Higgs particle and the 90 degrees n — n scattering cross section. We found that 
they are given by the product of a "universal" factor (identical for all regulariziations 
considered) which only depends on g and the dimensionless external momenta, and a 
non-universal factor that depends on the parameter bi but not on g or the momenta. 
The universal factor associated with the width is different from the one associated with 
the cross section, but the non-universal factors are identical for a fixed regularization 
scheme. These properties suggest that, to leading order in 1/N, all cutoff effects in on- 
shell dimensionless physical quantities (viewed as functions of dimensionless momenta) are 
given by an effective renormalized action, 

Seff = SR + cexp[ ]0 , (4.4.13) 

J J g 

where O is a renormalized operator, c is a ^ independent free parameter containing all 
the non-universal information (i.e. a function of bi only), and Sji is describing the usual 
universal part of physical observables with the unit of energy set by /tt = 1. In the 
general RG framework this representation of S^ff is not unreasonable if one accepts that 
at = cxD the number of independent operators that contribute to observables at order 
1/ decreases by one relatively to < 00. Because of the nonlinear relationship between 
c and the parameters in the bare action, the range in which c is allowed to vary depends 
on the type of action chosen. Different actions that realize the same c are indistinguishable 
to order 1/A^; all actions have large regions of total overlap but the triviality bounds are 
obtained from the area near the edges of the ranges in which c varies and therefore are 
somewhat dependent on the particular regularization scheme. Still, this dependence turns 
out to be weaker than what one would have been inclined to believe on the basis of the 
few simulations that had no continuous tuning ability for c built in. Thus, some of the 
stronger "universality" evident in (4.4.13) also makes its way into the triviality bound. 

We developed an approximate physical picture associated with the models we studied. 
When the regularized model is nonlinear one has to think about the Higgs resonance as 
a loose bound state of two pions in an / = 0, J = state. Pions in such a state attract 
because superposing the field configurations corresponding to individual pions makes the 
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state look more like the vacuum and hence lowers the energy. The four derivative term 
in the action can add or subtract to this attraction. We found that the smallest cutoff 
effects are obtained when the coupling bi of the four derivative term is set so that the term 
induces the maximal possible repulsion between the pions, postponing the appearance of 
the Higgs resonance to higher energies. 

Our findings concerning the Higgs mass triviality bound are summarized by figures 
7.3a - 7.3d and 7.4a - 7.4d for the four different regularization schemes considered. There 
the leading cutoff effects in the width and cross section are plotted versus the Higgs mass in 
TeV with the large N results presented for = 4. We decided to extract the bound from 
figures 7.4a - 7.4d since the cross section cutoff effects provide a more stringent criterion. 
For all regularization schemes we extracted the bound by restricting the allowed cutoff 
effects in the cross section to 4%, at center of mass energies up to twice the Higgs mass. 

To make predictions for the physically relevant case AT = 4 we used the known differ- 
ences between the N = oo and the N = A numerical results, for the simplest lattice actions, 
to extrapolate to the actions with four derivative terms that have not yet been studied 
numerically. We argued that the expansion in 1/N is expected to be "well behaved" in 
the region where the triviality bound is obtained, and because of that this extrapolation 
is sensible. 

Based on what we have learned, it seems that a more realistic and not overly conser- 
vative estimate for the Higgs mass triviality bound is 0.750 TeV, and not 0.650 TeV as it 
is sometimes stated. At 0.750 TeV the Higgs particle is expected to have a width of about 
0.290 TeV and is therefore quite strongly interacting. 
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APPENDIX A. 

In this appendix we show one example of the calculations needed to derive the asymp- 
totic expansions of the "bubble" integrals in (4.1). We shall work out explicitly only the 
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simplest bubble, /o,o(P^)- The other bubbles are calculated by similar means. Our problem 
is to calculate the first two terms in the expansion in of 

hfi{p ) - K{\p-k)K{\p + k) 

1 1 1 

4--E^^ (A.l) 



K{q) 0-2(2 + g2n) ^,2 „ Z-^ ^2 _^ 



TT 

Wj = exp[i(f>j] , (f>j = —{2j — 1 — n) for j = 1, 2, . . . , n . 



We rewrite the inverse propagator as 



1 " a 



— with WQ = . (A.2) 

j=0 



The values of the constants Cj can be read off (A.l). We now analytically continue in 
all the Wj, j = l,...,n, moving them to high positive imaginary parts. The bubble 
integral is now evaluated with this new propagator. We also analytically continue in the 
dimensionality d and first compute /o,o(P^) with the new propagator in d dimensions and 
then take the limit d ^ 4: 



hm y C,C. I 



d^ 



B{p) = hm V C,C. / . . , . (A.3) 



At the end we intend to analytically continue the tUj's back to their values on the unit 
circle. 

We intend to use Feynman's formula 

1 ds 
AB Jo [sA+{l-s)B]'^ ■ ^ ' ^ 

This is allowed as long as the straight line connecting the relevant points A and B never 
passes through the origin. This is ensured by all the tuj's being in the upper half of 
the complex plane and thus making the segments connecting A and B also stay in the 
upper half plane for all real values of the momenta. We introduce (A. 4) into (A.3) and 
do the angular part of the k integration. After a subsequent partial integration over the 
magnitude of k we obtain 



B{p) ^ lun C^Cj 



2ci7r5r(^) 0<f7<n * ^-^0 -^0 k'^ + s{l-s)p'^ + swi + {l-s)wj 

(A.5) 
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Let $s be the phase of Zs = s(l — s)p^ + swi + (1 — s)wj. It is clear that < $s < tt. 
We can deform the k integral from the positive axis to a ray starting from the origin at 
an angle This ray will always be in the first quadrant of the complex plane. The 

contribution from the arc at infinity is seen to vanish. Along the ray we can easily do the 
integral over k and obtain 



B{p)= lim ^ V C,Cj f\s[s{l-s)p'^ + swi + {l-s)wjt^ . (A.6) 



The phase of the integrand is ^^$5 so that the associated sheet of the logarithm is the 
first sheet cut along the negative real axis in the usual way. 

We now take the limit on the dimension and obtain, with an unambiguously defined 
logarithm, 



n 

2 



^ CiCj / ds / ds\og[s(l - s)p^ + swi + (1 - s)wj] 

'TT^ n/tf/„ ' Jo Jo 



(A.7) 



The first term vanishes because YlJ=o — 0- We separate out the terms containing wq 



B{p) = Ji + J2 + J3 
1 



Jl = 



1677^ Jo 

n „i 



/■I 2 

/ ds log[s(l — s)p 

Jo 



1 /"^ 
STT^n ^ Jo 

1 " rl 

J3= - ^Q^2^2 / log[s(l - + + (1 - • 



Since none of the integrals goes through the cut along the negative real axis they can be 
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expanded in p quite easily leading to 



Jl = 



167r2 



logp + 



87r2 



•^2 = ^E 



n r 9 



+ log - 1 



^3 



167r2 



[-2^2 



n r 9 



+ log tUj 



+ 0(^2)2) 

1 



167r2n2 



E 



WilogWi-WjlogWj 2/1 + 
^ - 1 +p I 



Wi — Wa 



+ 0((p2)2) 



WiWj 



2 - Wj)'^ {wi - Wj)^ 



(A.9) 



Finally we take the tuj's back to their original values, always staying on the first sheet, 
and never crossing the cut. We know then that in the above equations logtUj = for 
j >1 with the angles given by (A.l). The summations can now be performed by using 



1 77 — 1 

Z+WA ~ 1 + z^ 



A. 

dz 



log[l + z""] 



(A.IO) 



and derivatives of this identity. 



The calculation of the other bubble integrals to the same order needs some additional 
manipulations but is essentially of the same type as the one outlined above. 



APPENDIX B. 

In this appendix we shall evaluate the scaling violations in /5-functions defined along 
particular lines in the space of bare couplings of the simplest model 



(B.l) 



We shall work in the symmetric phase for simplicity. 



Introducing the aiixiliary field a — ^(j? we get a new action 



Si = -Trlog[K + ml + a] - / a 



(B.2) 
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Defining the bubble integral 



J(m2) 



1 



one obtains for the renormalized coupling 



(B.3) 



r;2 

9q 



(B.4) 



where the prime denotes diff'erentiation. The mass m? is given by the saddle point equation 

= ml+glj{m^) . (B.5) 

Suppose now that ^q(s), TOq(s) are a parametric description of a path in parameter 
space and we calculate along this path the variation of g with respect to m . We obtain 



dlogg^ 
dmn? 



= g^J"{m') + 



dloggQ 
drriQ 

2 T/f 



X 



X 



1 + g^J'{m^) 



(B.6) 



l + g^J'{m?) + J{m?)g^'^ 



It must be true that ^l^*^!" 7^ 00 because one cannot get to the critical surface without 

varying mg. Let us now consider the non-linear limit g^ ^ 00 mg — > — 00 with gf^jm^ — > ^. 
One can still have a path in parameter space ending at a critical point by varying ^ and it 
makes sense to compute j^f along this path. At ^fg = 00 one gets g^J'{vr?^ = —1 (see 
(B.4)) and therefore X = leading to 

dlogg'^ 



2 t/// z\ 
g J [m ) 



dm? 



(B.7) 



This answer is independent of the sequence of s dependent paths we took the large gfg, — Wg 
limit on, as it should be. 



Let us now compute in ordinary perturbation theory to one loop order. One gets 

d\oggl 



dm'' 



dmQ 



g'J{m')- 



dmg 



(B.8) 



with an explicit dependence on the sequence of paths, except if one happened to choose 
^^^1° = 0, but this is not necessary as the following example shows explicitly: 



2 1 
mg = - 



2A 



9l 



A 



(B.9) 
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The nonlinear limit is attained when A — > oo. We see that, unless a = 0, we shall have 
an answer that is different from the exact answer in (B.7) and that higher orders in (p' 
that have been neglected depend on a in such a way as to conspire to add up to give an ex. 
independent answer. The problem can be easily traced now to the fact that the quantity 
X in (B.6), which vanishes in the nonlinear limit, does so only as a result of keeping all 
the powers in g'^j'{rrP') that come in from the denominator. Since J' {w?) oc logm^ wc 
see that the culprit for finding a spurious a dependence in the leading order expansion of 
an a independent quantity is the truncation in the number of loops. As pointed out in 
section (4.3) graphs with an infinite number of loops contribute even at leading order in 
the physical coupling constant to the coefficient of the leading correction in inverse powers 
of the cutoff. All we are stressing here is that the unsolved problem of summing up the 
logarithms in the subleading corrections in the inverse cutoff is a serious problem. Maybe 
the 1 /N expansion is the right place to tackle this issue. 



In section 2 the linear model was discussed for the Pauli-Villars regularization and it 
was shown that the largest Higgs mass is obtained in the nonlinear limit of the model, where 
the four-point coupling A is taken to infinity. The action was restricted for simplicity to a 
three dimensional parameter space instead of the four dimensional one we really need to be 
concerned with (see introduction). In this appendix we will consider an action with four 
free parameters and again show that the largest Higgs mass is obtained in the nonlinear 
limit. However, our analysis is somewhat less general than the one in section 2 because it 
only covers a region close to the critical line. In section 2 we could do a better analysis 
because, with the particular Pauli-Villars regularization employed there, we had a simple 
closed formula for the "bubble" diagram, while here we work on the lattice and closed 
formulae for the "bubbles" are not available. We will use the F4 lattice to regularize the 
model. 

The action of the linear model with four bare parameters can be written as [6] 



APPENDIX C. 



S= - 




<x,x'> 



X 



X 



(C.l) 




To make the action bilinear on the ^ fields we introduce (as in section 2.1) two auxiliary 
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fields oj and a by inserting 



n r r d.{x)e^-^^^l^i^y^^ = 1 (C.2) 

^ J—oo ^'^ J— ioo 

into the functional integral. The zero mode of the ^ field is separated as in (2.1.7) and 
the N — oo pion propagator can be immediately read off. Prom the pion propagator and 
using standard conventions (see (3.3.5) and (3.3.6)) the pion wave function renormalization 
constant is found to be 

= i- . (C.3) 

As in section 2.1, to get the saddle point equations we integrate the H and tt fields 
and take the A?" — > oo limit. The saddle point equations are 

6k ^ ^ ^ (C.4) 

(1 - 2A - 24k) + 2(A + 12r]f)as + 3wl = ^i^s' , 

with 

.,,= ^-4 (C.6) 

and with Ji the F4 lattice integral defined in (5.1.1). The symmetric phase is characterized 
by = 0, LJs' > and the broken phase by v'^ > 0, Ug/ = 0. The phase diagram is obtained 
with the method described in sections 3.1 and 5.1. It can be drawn on a two dimensional 
plane spanned by the parameters a and b, where 



A + 12?7/ 
~(6^ 



i2 



(C.6) 



, ^ 1-2A A + 12?7/ 
^ = '-^-^«l6^ 

with ro = ^i(O) (see (5.1.5) ). There is a second order line described by 

b = hrQ, a > -2/iro, with h = ^^^^^3 . (C.7) 

This second order line turns into a first order line at the point a = — 2/iro, b = Htq. This 
point is a tricritical point. We see that the generic phase structure we obtained for the 
non-linear actions can also be found in the linear case. 

The physically relevant part of the phase diagram is the region near the second order 
line in the broken phase. In the broken phase the first of the saddle point equations (C.4) 
gives 

« = ^ (C.8) 
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where we have used = ^tt/^ (equation (4.1.3) ). We use this equation to trade the 
parameter k for /^r. 

Following sections 2.2 and 2.3 we obtain the propagators in Fourier space in terms of 
a matrix M as in (2.3.6). M is given by 



( 9{V)MI) 



M = 



\ 



0-s 





'Nv 
N 





N 



(C.9) 



N[X + 12?7/ + 3?7(Js - ^r]ig{p)] ) 



with /(p) = —j^logp^ + ci — C2P^ + 0{p'^logp^) (see (5.2.8)), the "bubble" integral defined 
in (5.2.7) and g{p), the kinetic energy term on the F4 lattice given in (2.5.5). 

Close to the second order line and using an analysis similar to the one in section 4.1 
we obtain equation (4.1.2) with the non-universal constant C depending on A, rj, and rjf 



m% ^ C2(A, r], r]i) exp[-967r2/^i?] 



(C.IO) 



where 



C{\ V, V) = exp 



47r^rQ 



cr2(A + 12r]/ + 3r]ac) 



(C.ll) 



In (C.ll) (7c = crc(A, r/, ry/), the critical value of as, is determined from equations (C.6,C.7) 
with Kc 



5av ■ 



From (C.IO, C.ll) it is clear that for a given mji the maximum gji = S-t^ is obtained 
when C(A, rj, rjf) is the smallest. If the factor A + 12rjl -\-2>rjac, appearing in the exponent in 
(C.ll), could take negative values then C would have been the smallest for those values. 
However, this factor cannot take negative values because of the second of eq. (C.7) which 
when written out is 



A + I2r]l + 3r7C7c > 



(C.12) 



This is a consequence of the fact that the physical region is the one close to the second 
order line and does not include the tricritical point or the region close to the "would be" 
second order line {h = hvQ, a < — 2/j,ro) which corresponds to an unstable saddle point. 
Therefore for given 77,77/ the minimum value of C(A, 77,77/) is obtained for A 00. We 
see, therefore, that the maximal Higgs mass is again obtained in the non-linear limit, at 
least as long as we restrict the observable cutoff effects sufficiently to justify the neglect of 
subleading cutoff effects in (C.IO). 



76 



REFERENCES 



U. M. Heller, H. Neuberger and P. Vranas Phys. Lett. B283 (1992) 335. 

U. M. Heller, M. Klomfass, H. Neuberger and P. Vranas, Nucl. Phys. B (Proc. 
Suppl) 26 (1992) 522; and work in progress. 

H. Neuberger, talk given at the workshop "Non-perturbative Aspects of Chiral Gauge 
Theories", Rome, 1992, to appear in the proceedings. 

H. Neuberger, in "Lattice Higgs Workshop", eds. Berg et. al.. World Scientific, (1988). 
H. Neuberger, Nucl. Phys. B (Proc. Suppl) 17 (1990) 17. 

G. Bhanot, K. Bitar, U. M. Heller and H. Neuberger, Nucl. Phys. B343 (1990) 467. 

H. Neuberger, Phys. Lett. B199 (1987) 536. 
S. Weinberg, Physica, 96 A (1979) 327. 
M. Einhorn, Nucl. Phys. B246 (1984) 75. 

M. Einhorn, in "Lattice Higgs Workshop", eds. Berg et. al.. World Scientific, (1988). 

C. Whitmer, "Monte Carlo Methods in Quantum Field Theory", unpublished Prince- 
ton University dissertation (1984). 

M. M. Tsypin, in "Lattice Higgs Workshop", eds. Berg et. al.. World Scientific, (1988). 
J. Kuti, L. Lin and Y. Shen, Phys. Rev. Lett. 61 (1988) 678. 
S. Coleman, R. Jackiw, H. D. Politzer, Phys. Rev. DIO (1974) 2491. 
H. Georgi, "Weak Interactions and Modern Particle Theory", Addison- Wesley, (1984). 
R. S. Chivukula, M. J. Dugan, M. Golden, preprint HUTP-92/A025. 
W. Bardeen, M. Moshe and M. Bander, Phys. Rev. Lett. 52 (,1983) 1188. 
M. C. Berger and F. David, Phys. Lett. B135 (1984) 412. 
M. Liischer and P. Weisz, Nucl. Phys. B318 (1989) 705. 

L. Lin, J. Kuti, Y. Shen, in "Lattice Higgs Workshop", eds. Berg et. al.. World 
Scientific, (1988). 

21) H. Neuberger, Phys. Rev. Lett. 60 (1988) 889. 

77 



22) A. Hasenfratz, K. Jansen, J. Jersak, C. B. Lang, H. Leutwyler and T. Neuhaus, Z. 
Phys.C 46 (1990) 257. 

23) R. Dashen and H. Neuberger, Phys. Rev. Lett. 50 (1983) 1897. 

24) H. Flyvbjerg, F. Larssen, C Kristjansen, Nucl. Phys. B (Proc. Suppl) 20 (1991) 
44; C. Kristjansen, H. Flyvbjerg, preprint, NBI-HE-90-55. 

25) M. Liischer and P. Weisz, Nucl. Phys. B290[FS20] (1987) 25. 

26) M. Liischer and P. Weisz, Nucl Phys. B295[FS21] (1988) 65. 

27) G. Bhanot, K. Bitar, U. M. Heller and H. Neuberger, Nucl. Phys. B353 (1991) 551. 

28) M. Gockeler, K. Jansen, T. Neuhaus, Phys. Lett. B273 (1991) 450. 

29) M. Gockeler, H. A. Kastrup, T. Neuhaus and F. Zimmermann, preprint HLRZ 92- 
35/PITHA 92-21. 

30) M. Gockeler, H. A. Kastrup, T. Neuhaus and F. Zimmermann, Nucl. Phys. B (Proc. 
Suppl) 26 (1992) 516. 



78 



